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Reduction of Oxides of Nitrogen 


in Vent Gases’ 


H. R. L. STREIGHT? 


Oxides of nitrogen are present in the effluent gases from nitric acid 
and nylon intermediates plants and if released in sufficient quantity they 
could be objectionable owing to their intensive reddish-brown color and 


to their corrosive and poisonous nature. 


A paper study showed that the 


oxides of nitrogen could be reduced by absorption in an alkaline solution. 
Based on pilot plant tests, a large fume abatement unit was designed, 
installed, and operated in a satisfactory manner to reduce the concen- 


tration of the oxides of nitrogen. 


In 1956, the color of the combined 


effluent gases was practically eliminated when a more economic waste 
ammonia-caustic liquor replaced the dilute caustic scrubbing solution 
and finely dispersed ammonium salts were removed in a Venturi scrub- 
ber. The method may be used to eliminate the objectionable brown- 


ish-red color of effluent gases leaving a nitric acid plant. 


D Pont of Canada operates a 
complex organic chemical plant 
at Maitland, Ont., for the production 
of nylon intermediates, adipic acid 
and hex amethylenediamine, the raw 
materials which are combined and 
polymerized to make nylon filament 
and staple fibre at its Kingston plant. 
When it was decided to erect the 
Maitland plant which is located in 
a recreational and agricultural area, 
an extensive pollution abatement pro- 
gram was started and continued after 
the plant was in operation, to deter- 
mine the amount of certain contami- 
nants in the effluent gases and to 
learn if they could be reduced be- 
fore release to the atmosphere. 

The stack gases from a nitric acid 
plant, using the ammonia oxidation 
process, are released as brownish-red 
plumes as are reaction off-gases and 
vent gases from tanks and equipment 
in the nylon intermediates plant. The 
objectionable color is caused by a 
concentration of nitrogen dioxide 
(NO,) above 50 p.p.m. On oxida- 
tion, nitric oxide (NO) forms objec- 
tionable NO,. Nitrogen dioxide can 


react with water vapor or rain to 


form nitric acid. According to re- 
cent studies on pollution, the oxides 
of nitrogen can affect atmospheric 
pollution in an indirect way since 
they act as catalysts and oxidize or- 
ganic materials in the presence of 
ultra-violet light or bright sunlight 
to objectionable compounds. The 
American Petroleum Institute and 
other scientific organizations in Los 
Angeles are investigating this prob- 
lem. The other constituents of the 
stack gases are nitrous oxide (N,O), 
carbon dioxide, oxygen and nitrogen. 
The N,O does not cause any anxiety 
as it is a colourless gas occurring in 
such dilute concentration that its 
anaesthetic properties would be un- 
noticed. Oxides of nitrogen hence- 
forth refer to NO and NOs. 

The recovery of low concentra- 
tions of NO and NO, is not 
economic as the gases are too dilute 
for conversion into nitric acid. The 
problem accordingly was to remove 
a potential nuisance by reducing the 
emitted oxides of nitrogen to such 


a value that after dispersion from a 
stack, their maximum ground con- 
centration would be less than the 


1Manuscript received November 8, 1957. 


2Principal Chemical Engineer, Engineering Department, Du Pont Company of Canada (1956) Limited. 
This article is based on a paper presented at the 40th Annual Conference of The Chemical Institute 


of Canada, Vancouver, B.C., June 3-5, 1957. 
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safe threshold limits accepted for in- 
dividuals and vegetation; then there 
would be no injury to the health of 
persons or damage to the apple trees 
in adjacent orchards and flowers in 
nearby gardens. Since the maximum 
safe ground concentration of each 
contaminant varies with the individ- 
ual and the actual plant, the follow- 
ing values were accepted as the safe 


and maximum ground _ concentra- 
tions. 
TABLE 1 
THRESHOLD LIMITS FOR OXIDES OF 
NITROGEN 
Max. Ground Level 
Concentration 
NOs, p.p.m. Duration 
Persons 5 (1) continuous 
10 (3) 8 hours/day, 
5 days/week 
Vegetation <25 (2) continuous 


According to Table 1, the concen- 
tration of the NO and NO, released 
from a stack should be such that 
after dispersion in the atmosphere, 
the maximum ground concentration 
should not exceed 10 p.p.m. At this 
concentration, the effluent gas would 
be colorless as the reddish-brown 
color is clearly visible only above 


50 p.p.m. NO,,. 


Literature search 

The technical literature was re- 
viewed to provide information on 
the properties and reactions of NO 
and NO,, and methods for their 
reduction and removal. 

Nitric oxide is a colorless and re- 
latively stable gas which is only 
sparingly soluble in water. It slow ly 
combines with oxygen to form the 
reddish-brown nitrogen dioxide 


3 





which reacts with water according 
to the following equation: 


3NO) + H.0 — 2HNO; + NO 


When NO and NO, occur together, 
they establish an equilibrium with 


the trioxide (N, O,); NO, is also in 
equilibrium = w ith ‘the tetraoxide 
(N.O,) 


a constants between the 
various oxides of nitrogen are shown 
in Table 2. 


TABLE 2 
EQUILIBRIUM CONSTANTS 


Equilibrium 


Reaction Constant 


K, at 20°C. 
2NO>2 >NO; 10.1 
2NO +0, = 2NO, 4.2x 10'? 
The above equilibrium constants 


were applied to a gas containing 
0.05% NO and 0.025 volume % 
NO.,,. Quantitative calculations show- 
ed that the concentrations of tri- and 
tetraoxides would be exceedingly 
small and could be neglected. As the 
gas would contain approximately 4% 
oxygen, the third reaction would 
give substantially complete yields of 
NO, provided sufficient time was 
available; the effluent gas would re- 
quire more than 30 minutes to re- 
duce the initial NO content by 50% 
and more than 5 hours for 90% 
Hence oxidation would be too slow 
for practical purposes, even with 
additional oxygen. A _ satisfactory 
catalyst for this reaction has not 
been found although long sought by 
nitric acid manufacturers who are 
always interested in improving the 
efficiency of their absorption units. 

The following methods of  re- 
ducing NO and NO. have been 
studied, although often not exhaus- 
tively. 

1. Absorption (Absorbers, 
lation Columns and 
Scrubbers ) 

(a) water 

(b) alkaline solution 
(c) alkaline sulphite solution 
Absorption 
Combustion 
Chemical reaction, 
following materials: 
(a) sulphur dioxide 
(b) chlorine 


Distil- 
Venturi 


> w nN 


with one of 


1. Absorption. For cheapness and 
availability, water was first con- 
sidered as the absorbing medium in 
absorbers and Venturi — 
The absorption of NO and NO, 

water has the objection of aan 
only two-thirds of the NO,. The 


fd 


balance is converted to NO which, 
being insoluble in water, is released 
to the atmosphere where it oxidizes 
to NO,, forming the objectionable 
reddish-brown color. 

Alkaline absorption has been used 
to reduce the last traces of oxides 
of nitrogen leaving a Norwegian 
nitric acid plant, and to clean the 
fumes from an electroplating piant 
and recently from a Canadian chemi- 
cal plant“); unfortunately information 
was not recorded for the design of 
another plant. The rates of absorp- 
tion of NO and NO, in alkaline 
solutions have been determined in 
several laboratories but not in suffi- 
cient detail for design. The mechan- 
ism of absorption of low concentra- 
tions of NO, in caustic solutions was 
studied by Chambers and Sherwood) 
and their results confirmed by E agle- 
ton, Langer and Pigford). The rate 
of absorption was found to be gas- 
film controlled and the absorption 
coefficients were considerably lower 
than expected. Various mechanisms 
were proposed to account for this 
discrepancy. The absorption rate in 
caustic soda was less than in water; 
the coefficients varied with the 
liquor rate, gas rate and water vapor 
content of the gas. 

Further experiments using N,O,, 
instead of NO,, showed that N. O, 
reacts rapidly with NaOH probably 
according to the following reaction: 


N2O; + 2NaQH — 2NaNO, + H20 


The formation of N,O,, from NO 
and NO., in a distillation column 
was shown by the increase of plate 
efficiency in proceeding up the 
column. Tests were not made with 
an NO and NO, feed stream™, If it 
is assumed that absorption proceeds 
through the formation of N,O,, the 
molal concentration of NO must not 
be greater than the molal concentra- 
tion of NO, — otherwise excess NO 
will be released to the atmosphere 

After our test work was com- 
pleted, other experimenters showed 
that removal was improved by sup- 
plying a maximum gas-liquid contact 
area applicable to absorbers, distilla- 
tion columns and Venturi scrubbers 
(water is sprayed axially into the 
gas stream passing at a high velocity 
through a Venturi throat). For low 
concentrations of oxides of nitrogen 
(1% or less) spray towers and 
bubble cap columns were found to 
have low removal efficiencies when 


water was used as the scrubbing 
medium (°: 9), 
Little information was found on 


scrubbing with an ammonium hy- 
droxide solution although mention 
was made of the formation of am- 


nitrate and ammonium ni- 
trite in the scrubbing liquor and of 
preferential reaction of the liquor 
with NO, over NO. A strong alka- 
line solution of a common sulphite 
would not be satisfactory for scrub- 


monium 


bing the Maitland gases as the 
oxygen would soon oxidize the sul- 
phites to sulphates. 

2. Adsorption. The adsorption of 
oxides of nitrogen by silica gel (or 
other adsorbents) is satisfactory for 
NO, concentrations — higher than 
0.1% and for NO concentrations 
higher than 1 to 1.5% (11). Silica 
gel containing adsorbed NO, cata- 
lyzes the oxidation of NO and the 


formed NO, is removed by adsorp- 
tion on the silica gel w hich on heat- 
ing produces a higher concentration 


of NO, that can be used in a nitric 
acid plant. The equipment is expen- 
sive and operation costs would be 


high for low concentrations of the 
oxides of nitrogen. 


3. Combustion. The effluent gases 
could be treated with a small amount 
of hy drogen over a suitable pallad- 
ium catalyst. The presence of 
oxygen in the Maitland effluent gas 
obviates this method of removal 
since, in all likelihood, the hydrogen 
might react with the oxygen before 
the oxides of nitrogen (12, 18), 


4. Chemical Reaction. The oxides 
of nitrogen can combine with sulphur 
dioxide or chlorine. Chemical reac- 
tion does not appear attractive as 
sulphur dioxide would have to be 
added to the effluent gases at a high 
temperature and the rate of reaction 
is very slow either for sulphur di- 
oxide or chlorine (14, 15, 16), 
Proposed Method. Based on _ the 
literature search, the most promising 
method appeared to be absorption 
in a sodium carbonate or hydroxide 
solution. There was a potential ob- 
jection owing to the presence and 
concentration of CO, in the com- 
bined effluent gases. Carbonate and 
bicarbonate would be formed in the 
solution until a steady state was 
reached when the carbonate/bicar- 
bonate ratio and the total ionic 
strength would generate a back pres- 
sure of CO, comparable with the 
forward pressure in the transpired 
gas“), Accordingly a_ pilot plant 
absorber system was recommended 
as the next step in determining the 
potentially useful concentrations of 
alkaline scrubbing solutions and ab- 
sorption efficiencies, for low concen- 
trations of NO and NO,. The test 
runs would use sodium carbonate 
solutions after establishing the 
method with water. 
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A Venturi scrubber might have ‘a 
lowe: capital cost than an absorber 
for removing oxides of nitrogen. Ac- 
cordi ngly pilot plant tests were 
scheduled on a Venturi scrubber to 
determine the number that would 
be required. 


Pilot plant 

An absorption column was con- 
structed from 12-in diameter sections 
of glazed pipe containing a 5 or 10- 
ft. depth of random packed Raschig 


rings. The characteristics of the 
packing were as follows: 
Free space —77% 
Ring dimensions —1 x 1x 0.1-in. 
Liquid feed rate, 
for complete 
wetting —1.0 cu. ft./min. 
Each inlet gas line was equipped 
with a rotameter, manometer and 
facilities for gas sampling. Dilute 


air/NO, mixtures were made to the 
correct concentrations. Liquid NO, 
was evaporated by an immersion 


heater oper rated at controlled wat- 


tage inside a Dewar container. This 
evaporator was calibrated empiri- 
cally and evaporated NO, was 


measured before being mixed with 
the main air stream. A pump recir- 
culated the liquor through the ab- 
sorber, countercurrent to the gas 
stream. Special sampling and analy ti- 
cal techniques were required largely 


because of the low levels of con- 
centration of these components in 
the waste gas. The collection of a 


large volume for sampling, coupled 
with a combination of photoelectric 
and chemical analytical methods, 
proved quite satisfactory. The an- 
alytical methods were revised and 
improved throughout the  experi- 
mental work. 

The initial experimental work was 
conducted in the period April-Sep- 
tember 1952, and involved the pre- 
paration of a “synthetic” waste gas 
containing 0.05 volume % NO, 
which was absorbed by water in the 
pilot plant absorber. The results of 
this first group of runs showed that 
the absorption of NO, was limited 
by the gas-film coefficient and the 
“Height of a Transfer Unit” 
(H.T.U.) for absorption of NO, i 
water was much higher than pre- 
dicted on the basis of the two-film 
diffusion mechanism, by a ratio ap- 
proximately 33:1. Several runs were 
conducted with NH, and Cl, present 
as contaminants in the waste gas, in 
concentrations of 0.05-0.10 volume 
%. The H.T.U. values obtained with 
NH, were comparable to those with 
uncontaminated waste gas, while 
those obtained with Cl, were even 
higher. 
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Subsequent runs centered around 
the absorption of NO, in a solution 
of sodium carbonate. Complications 
were expected with the use of this 
solution as the waste gas contained 
approximately 2.4 volume % of CO,,. 
Recirculation of the absorbing | iquor 
resulted in a continuously increasing 
bicarbonate/carbonate ratio in the 
liquor (up to the point of equili- 
brium), with an attendant decrease 
in the absorpiton efficiency of the 
absorber. The absorption efficiency 
of NO, in a gas containing 2.4 

CO, dropped from 40- 


volume % 

50% for a bic arbonate/carbonate 
ratio of 0-0.3, to 20-25% for a bi- 
carbonate/carbonate ratio of 4. The 
possibility of ventilating the recircu- 
lated liquor in order to reduce the 
bicarbonate/carbonate ratio was in- 
vestigated in a series of runs. Al- 
though the _ bicarbonate/carbonate 
ratio did decrease, the efficiency of 


TABLE 


oxide absorption failed to show the 
expected recovery. Thus regenera- 
tion of the recirculating liquor in 
terms of CO, content, was ruled out. 
One important aspect in the utiliza- 
tion of the alkaline absorbent is the 
extent to which the absorbent can 
be exhausted, in terms of sodium 
ion content, before it becomes use- 
less. Exhaustion is defined as the 
percent of sodium ions originally 
present in the solution which have 
been converted to nitrate and nitrite 


salts. Initial work was conducted 
with 50% exhausted liquor contain- 
ing equilibrium amounts of bicar- 


bonate and carbonate, as well as 50% 
exhausted liquor containing higher 
than equilibrium _ bicarbonate/car- 
bonate ratios. 

A comparison of water and 


Na,CO, solutions as absorbing media 
is shown in Table 3. 
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COMPARISON OF WATER AND NasCO; So_utions As NO» ABSORBING MEDIA 


Gas Rate 60 cu.ft./min.* 


Water 


Pco, O atm. 


Pyo, 
21° 


Absorption efficiency....... ‘ 


*Tests were made on gas rates ranging from 40 to 100 cu.ft. 


the calculated carry-over rate. 


(93% 


Packed Depth = 5-ft. 


0.0005 atm. 


Liquor Rate 1 cu.ft./min. 


57 N las 2COs Solution 
0. 024 atm. 


0.0005 atm. 


at zero exhaustion = 50% 
at 50% exhaustion = 22% 
‘min., 118 cu.ft./min. being 


EXHAUSTED AT START) 


(85% EXHAUSTED AT START) 
NO, 


ONLY 


' 
90 


? 
92 


~ EXHAUSTION 





Figure 1—Effect of exhaustion on H.T.U. 
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The tests showed that a 50% ex- 
hausted alkaline solution offered only 
slightly better absorption efficiency 
than water. The disadvantage of 
water was the very slight solubility 
of NO and its release to the atmos- 
phere. Thus as the waste fumes con- 
tained NO and the absorption of 
NO, in water would liberate one- 
third of its volume as NO, sodium 
carbonate solutions were selected as 
the absorbing medium. 

During the period December 1952- 
February 1953, the experimental data 
on carbonate solutions were ex- 
tended to higher exhaustion levels, 
using synthetic waste gases contain- 
ing NO.,,. The absorption coefficients 
for liquor exhaustion levels of 78- 
100% were somewhat scattered, pos- 
sibly because of impurities in the 
prepared “waste” gas. As Figure 1 
will illustrate, the absorption effi- 
ciency remains relatively constant 
up to about 90-95% liquor exhaus- 
tion. The H.T.U. values were cal- 
culated assuming zero back-pressure 
for the dissolved NO,,. 

Table 4 arbitrarily divides the tests 
of Figure 1 and shows the observable 
trend. 


TABLE 4 
THE EFFECT OF EXHAUSTION ON 
H.T.U. For NO» 
Exhaustion | H.T.U. | Devia- | No. of 
tion Results 
78— 88% 16.5 ft. +4.8 ft 9 
88— 92% | 17.4 Le 8 
92— 96% | 20.8 5.0 2 
25.4 9.4 10 


96-100% | 


Before designing a large absorber, 
it was important to learn if HM... 
was enhanced in an absorber with a 
short depth of packing, 5-ft., by ab- 
sorption at the top and bottom ends. 
The effect of packed depth on the 
absorption efficiency of alkaline solu- 
tions was examined by performing 
a number of runs with a 10-ft depth. 
The observed efficiencies were twice 
those obtained with 5-ft. of packing. 

The effect of pH on absorption 
efficiency was studied by using a 
dilute caustic solution as the absorb- 
ing liquor. The CO, in the waste 
gas rapidly converted the solution 
to one of sodium carbonate. The pH 
of the solution, initially 12.8, de- 
creased to a value slightly higher 
than 7, providing an opportunity for 
continuously testing absorption effi- 
ciencies in this pH range. No marked 
variation of efficiency was noted. 

At this stage of the experimental 
work, additional process information 
indicated that the concentration of 
NO in the combined effluent gases 
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TABLE 5 
NO ON THE 


THE EFFECT OF 


Run . 

No. Exhaustion I NO, 

52 50% 0.001 atm. 
53 50% 0.001 

59 55% 0.001 

60 §5% 0.001 

57 91% 0.001 

58 91% 0.001 


from the plant, previously assumed 
to be of a molal ratio of 2NO/1 
NO,, would more likely be 1/1 or 
even slightly lower in NO. This is a 
more favourable ratio for absorption 
in sodium. carbonate/ bicarbonate 
solutions which may take place ac- 
cording to the following reactions: 


NO + NO: @ N,O; 


N.0; + 2 NaHCO; > 2 NaNO» + 
2 CO. + H.O 


O,; a Na.CO; 2 NaNO» + CO, 


Initial work with mixed oxides of 
nitrogen in the ratio of 1 NO/1.1 
NO, showed an improvement in 
over. rall oxide absorption efficiency, 
in agreement with the findings of 
earlier workers ('%, 1"), Additional ab- 


sorption runs were made on syn- 
thetic “waste” gas containing mix- 
tures of NO and NO, in a 1/1 ratio, 


with absorbing liquor exhausted to 
50% and 90%. The H.T.U. values 
are based on total oxides of nitrogen, 
assuming no back pressure of dis- 
solved oxides. The H.T.U.’s were 
somewhat erratic as shown in Table 5. 
The only reason which can be ad- 
vanced for this inconsistency is the 
possibility that the oxides of nitro- 
gen, produced in the laboratory for 
runs 57, 58, 59 and 60, contained 
impurities. 

For 
solubility 
_ on the 


plant design a number of 
determinations were con- 
system NaHCO, - 

Na,CO, - NaNO, - NaNO, - H,O, 
the liquor samples being saturated 
with a gas comparable to the com- 
bined Maitland effluent. Experiments 
showed that at liquor exhaustions of 
over 30%, precipitation would not 
occur. In the range 0-30% exhaus- 
tion, NaHCO, precipitated at a tem- 
perature below 8.5 and 5°C. for 0 
and 20% exhausted solutions. The 
solutions easily became super-satur- 
ated and even after seeding with 
NaHCO,, there was a_ protracted 
time before crystal growth gave a 
visible precipitate. 


was made _ be- 
1952, and Febru- 


An investigation 
tween December, 


ABSORPTION OF NO» 


Standard No. of 
Pxo H.T..U. Deviation Results 
14.1 ft. +1.6 : 
0.001 | 0.9 6 
- 1.8 0.1 2 
0.001 V4 1.8 4 
16.9 eis 4 ™ 
0.001 20.6 5.3 7 


ary, 1953, on using a Venturi scrub- 
ber for absorbing the oxides of 
nitrogen in a 5% sodium carbonate 
solution. Liquor was introduced into 


the throat of a Venturi tube by 
radial and axial injection methods. 
Nitrogen dioxide alone, and a 1/1 


ratio of NO and NO.,, were passed 
through this scrubber at absorption 
efficiencies of less than 40%. As a 
consequence Venturi scrubbers were 
not considered practical for a plant 
installation. 


From the preceeding description 
of experimental work, it was con- 
cluded that the absorption of the 
oxides of nitrogen could be most 
effectively performed in an alkaline 
solution of sodium carbonate or hy- 
droxide. An absorption column de- 
signed for the removal of NO, will 
be equally capable of dealing with 
NO/NO,, mixtures up to a 1/1 mole 
ratio although the tests did not con- 


sistently reproduce the well-docu- 
mented claim that N,O, is more 
easily absorbed than NO, in alkaline 


media. For design, the H.T.U. can 
be safely taken as equal to 20-ft. for 
random 1-in. rings, up to 90-95% ex- 
haustion of alkalinity; the absorption 
efficiency will decrease if the liquor 
is exhausted more than 90-95% 


Fume absorption installation 


Basis of Design. In 1953 a fume 
absorption system was designed, 
the basis of the pilot plant runs, to 
reduce the oxides of nitrogen in the 
effluent gases of the nitric acid and 
nylon intermediates plants in a satis- 
factory manner for the forecasted 
and maximum loads. The forecasted 
load represented normal operation 
after plant start-up and the maximum 
rate corresponding to an increased 
production of nitric acid. 

Vent gases containing NO 
NO, come from three sources: 


and 


i Off-gases from the nitric acid 
plant. 


2. Reaction off-gases. 
Vent gases from tanks and pro- 
cess equipment. 
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Several schemes were investigated to 
determine whether all or part of the 
above effluent gases should be passed 
through one or more fume absorbers. 
Dilution by the addition of air to 
the vent fumes having the highest 
concentrations of NO and NO, was 
considered and rejected as pollution 
would be too high in the immediate 
vicinity under adverse meteorologi- 
cal conditions. In the adopted 
scheme, all efHuent gases containing 
NO and NO, were passed through 
a single absorber. 


The concentrations of NO and 
NO, were calculated on the best 


available information, attention being 
given to the increase of NO and 
NO, content of the off-gas with in- 
creased production and higher cool- 
ing water temperature. As the esti- 
mates made on the vent gases from 
tanks and other equipment vary 
widely with changes in operation, an 
allowance was made for surges. Since 
start-up, the flow and composition 
of the majority of the off-gases have 
altered, owing to the development of 
economic process changes. The pre- 
sent flows differ from the values used 
for design. 

The NO and NO, contents of the 
off-gases were calculated for various 
depths of packing in the fume ab- 
sorber. The proposed absorber with 
a 40-ft. depth of wetted packing, 
equivalent approximately to 1.2 
H.T.U., should give conservative 
exit concentrations of NO and NO.,, 
particularly since calculations show- 
ed that diffusion would keep the 
maximum ground concentration be- 
low the threshold limit. The results 
are incorporated in Table 6 and 
should be compared with the thres- 
hold limit of 10 p.p.m. 

The design flows and compositions 
of the inlet and off-gases are given 
in Table 7. 

The liquor system was designed 
for batch operation, starting with 
3.5-4.0% caustic solution which 
would be essentially exhausted be- 
fore the spent liquor was sent to 
sewer and replaced with fresh liquor. 


Equipment. The fume absorp- 
tion system was designed in 1953 to 
provide the greatest possible reduc- 
ticn of the oxides of nitrogen in a 
single absorber with ancillary equip- 
ment. This system is shown in 
Figure 2. 

The low pressure off-gases are col- 
lected in appropriate process build- 
ings and passed through a blower to 
the fume absorber. High pressure 
reaction off-gases, from a nylon in- 
termediates process building and the 
nitric acid plant, pass directly to the 
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TABLE 6 


NO anp NO» CONCENTRATIONS 


Gas to absorber. 


Gas from absorber. 


Max. ground concentration,* based on design 
of absorber and after diffusion in a 2 to 3 


mi./hr. wind. . 


*The maximum ground concentrations have been found to be 50 to 67% of the calculated 


values (20). 


TABLE 7 


Forecast Load 
p.p.m. 


5,200 
1,850 
| 


0.49 (av.) 
1.2 (max.) 


FLOWS AND COMPOSITIONS FOR FUME ABATEMENT SYSTEM 


Total Inlet Gas 


Forecast 


Load 
Flow, Ilb./hr. es : 19,417 
Flow, GuGb fT so aces oe dsw 247 ,600 
‘Temperature, °C... cies 43 


Pressure, Ib./sq.in.ga. 


Density, MR ACU cic cerca 0.078 

Composition, Ib./hr. ; 
Nitrogen is |} 17,074 
Oxygen. 658 
Water... ae 117 
Carbom dioxide. ...5..5 05.6056. 696 
N.C 771 
NO** 23 
DON owed, “er Cereals 78 


* Measured at operating conditions. 


Design Conditions 


Maximum Forecast 
Load Load 
31,931 19,572 

410,700 250,000 
43 25 
2 1.3 
0.078 0.078 
28 , 268 17,074 
1,054 658 
156 339 
870 694 
964 771 
139 8 
480 28 


€ 


** Variation per inlet effluent stream, 0.05 to 33 wt. %. 
***Variation per inlet effluent stream, 0.07 to 49 wt. %. 


TABLE 8 


Liquor SYSTEM 


Batch size, 
Circulation rate, 


Time/batch, hr. 

Forecast load. 

Maximum load... 
NaOH exhaustion, aces wtijecian 
Composition of purge liquor 

NaNO» wt. % 

NaNOs, wt % 
absorber in which the combined 


gases are scrubbed by a caustic soda 
so!ution before being vented to the 
atmosphere. The knockout drum for 
recovery of condensed nitric acid 
and the separator for removing 
vapor from the circulating liquor are 
positioned in the circulating system 
as shown in Figure 2; later equip- 
ment changes for use of the system 
with a waste ammonia-caustic soda 
solution are marked separately, 


RN csr e ta wees wh 
Ue Ba Pas cis cme Senn 
Fresh batch composition, wt. % NaOH.............. 


Rw 
wun 
woo 


Maximum Load 


1.9 (av.) 
4. 


7 (max.) 


Off-Gas 


Maximum 


Load 


Oo = 


— in 


WMrm bdo 


31 


2 


, 893 
412,000 
5 


1.3 


0.078 


28 
1 


Nw 


On 


, 268 
054 
557 
856 


In designing the absorber, the fol- 
lowing points were considered: 


(a) A suitable packing. 


(b) A high gas velocity to im- 
prove absorption efficiency. 
(c) A liquor rate to wet all of the 


packing. 


(d) H.T.U. based on the results 
of the pilot plant tests cor- 


rected for 


changes 


in 


packing and gas velocity. 


the 
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Figure 2—Fume abatement system. 


(e) A packed tower height com- 
mensurate with optimum in- 
cremental removal of NO and 
NO,,. 

The 8-ft. diameter absorber con- 
tains a 41.5-ft. depth of 2-in. random 
packed Raschig porcelain rings, 18- 
in. acting as an entrainment separa- 
tor around the 26 liquor dip- pipes. 
The top of each dip-pipe contains a 
sloping orifice disc designed for a 
head of 9-in. of liquor and the bot- 
tom is flared and placed over a splash 
block located in the rings. The inlet 
liquor pipe divides inside the ab- 
sorber and feeds a single central 
liquor feed trough w hich’ passes the 
liquor through notches into five 
liquor distributors, with attached dip 
pipes, supported across the absorber. 

The absorber was constructed of 
rubber lined mild steel. Transfer 
lines and all other equipment which 


8 


come in contact with the waste gas 
were constructed of stainless steel 
(type 304). The equipment and lines 
handling the alkaline solution were 
constructed of mild steel. 

The packing was dumped into a 
water filled section of the absorber 
and can be removed through a man- 
hole above the supporting grid hold- 
ing the packing. The distribution of 
liquor, with no gas flowing through 
the absorber, can be seen through an 
open manhole under the grid. The 
fumes enter near the base and leave 
at a high velocity through a 15-in. 
diameter stack attached to the top 
centre of the absorber. 

A number of safety precautions 
are incorporated in the design. For 
instance, all headers leading from the 
tanks to the waste gas blower, are 
supplied with vents which are open 
to the atmosphere. These vents have 


the dual function of allowing the 
escape of the waste gas to the atimos- 
phere, in the event of a blower stop- 
page, and preventing the blower 
from pulling a vacuum on the tanks 
by admitting atmospheric air to the 
headers in the event of a waste gas 
flow decrease. Broken rings are re- 
tained by wire mesh strainers on the 
pump suction lines and stopped by 
an internal wire covered perforated 
pipe through which the liquor passed 
before leaving the circulating tank. 
No pump trouble has been experi- 
enced although many finely-broken 
rings were removed in the first few 
weeks of operation. On failure of 
a circulating pump or during a 
dump, a mercoid switch stops the 
blower and opens an_ automatic 
valve, venting the high pressure re- 
action off-gases to a stack and pre- 
venting the fumes from entering the 
absorber. At the same time an alarm 
on the blower sounds and the opera- 
tor vents the waste gases from the 
nitric acid plant through another 
stack to prevent acidic conditions in 
the absorber. 

Plant Tests. Although a caustic 
soda liquor proved satisfactory in 
reducing NO and NO, in the vent 
gases, the operating costs could be 
reduced substantially by using an 
available waste ammonia-caustic soda 
solution (0.5% free NaOH, 0.45% 
Na,CO, and 0.03-0.5% NH,) which 
was tested during 1955 in the fume 
abatement system as an alternative 
scrubbing medium. The waste stream 
was cooled to 35°C. to avoid damag- 
ing the rubber lining of the absorber 
and passed through ‘additional piping 
to the fume abatement system. Several 
plant tests of a semi-continuous 
nature were made to determine the 
percentage exhaustion, and the color 
and appearance of the plume leaving 
the absorber. The utility of the am- 
monia-caustic liquor was better than 
80%. An analysis of the solution 
showed a straight line relationship 
between pH and the degree of ex- 
haustion of the spent liquor between 
40 and 100% exhaustion. On 100% 
exhaustion, equivalent to a pH of 
7.5, the liquor released CO, gas. This 
gave a ready method of control since 


a pH _ meter, with recorder and 
alarm, could replace an operator 
who formerly spent considerable 


time on control analysis. Based on 
preliminary analytical tests, the re- 
moval efficiency of NO and NO, 
with the waste liquor was as good as 
that obtained with a caustic soda 
solution. This comparison is not out 
of line with the early pilot plant tests 
conducted with NH, present in the 
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gas phase. The plume from the stack 
of the fume absorber changed from 
a light reddish-brown to an intense 
white color; it is quite possible that 
the whiteness of the plume obscured 
any brown coloration of the gas. 
The appearance of the white color 
was attributed to the presence of 
aminonium nitrate and nitrite re- 
leased in a finely divided state; the 
nuisance that can be caused by these 
white solids has been reviewed in a 
recent article ?"), This objection had 
to be overcome before the waste 
liquor could be used for continuous 
operation. 

The addition of NH, to the base 
of the absorber and subsequent in- 
stallation of water sprays above the 
packing failed to reduce the white 
plume which under certain meteoro- 
logical conditions tended to settle 
on the ground downwind from the 
stacks. 

A microscopic examination of the 
solids showed that they were very 
small, 0.1 to 10 microns in size. As 
wet scrubbers are satisfactory for 
removing such small particles, tests 
were made in a pilot plant Venturi 
scrubber which was reported to give 
satisfactory removal for sub-micron 
sized particles. Information on the 
design and applications of the Ven- 
turi scrubber is given by Johnston 
and Roberts (22). As the inlet water 
enters the throat of the Venturi and 
is dispersed into small particles by 
the high velocity of the gas, intimate 
mixing is maintained between the 
liquid and the gaseous phases. When 
water was added to the test Venturi 
scrubber, 60 to 90% of the white 
fumes was removed, higher values 
being obtained when steam was 
added to coalesce the particles prior 
to the scrubber. 


System improvements. Following 
the test work conducted with a Ven- 
turi scrubber for the removal of fine 
solids, the findings were reviewed 
with Sheldons Engineering Limited, 
Canadian suppliers of Venturi scrub- 
bers. Their proposal for a Venturi 
scrubber and separator was accepted. 
Figure 3 illustrates the Venturi 
scrubber and separator installed in 
the fume abatement system. A simple 
separator follows the scrubber to re- 
move large liquid droplets, 50 to 150 
microns in diameter. 

The Venturi scrubber and separa- 
tor were designed for a gas flow of 
7,000 cu. ft./min. and a liquid flow 
of 42 U.S. gal./min. with a total 
pressure drop of 22 in. W.G. and 
recovery efficiency of 99.0% plus. 
The dimensions of the Venturi 
scrubber are 9-1/4 and 18-in. for 


throat and end diameters and 10-ft. 


6-in. for length and of the separator, 


4-ft. for diameter and 9-ft. for length. 


Performance 


The fume absorber’ was first 
checked for liquor distribution un- 
der varying rates of flow. Water 
was found to leave the packing in 
a uniform manner across the tower, 
without any sign of cascading down 
the wall. Good distribution was 
again noted after a year’s operation 
although an additional 6-in. of pack- 
ing had to be added to retain its 
initial depth. 

After the start-up of the plant, a 
feed liquor containing 3.5-5% caustic 
was circulated in the absorber until 
near exhaustion when it was replaced 
with a fresh batch. Because of an- 
alytical difficulties no reliable ab- 
sorption data are available for this 
period of operation, although the 
color of the off-gas stream indicated 


substantial reduction of the NO, 
content. 

Adequate material balances have 
been developed around the absorber 
with the waste ammonia-caustic soda 
solution as the scrubbing liquor. 
These studies were conducted at a 
time when the plant was operating 
at a capacity halfway between the 
forecast and maximum loads. The 
quantities of NO and NO, present in 
the waste gases were considerably 
reduced by improvements and 
changes in plant operation. In Table 
9 the plant results are compared with 
design figures derived from Table 7. 

Owing to surges, the quantity of 
oxides of nitrogen in the waste gas 
varies considerably. 

The caustic scrubbing solution was 
circulated on a _ batch basis with 
periodic titrations and dumped when 
it had attained a 95% exhaustion. - 
Since employment of the waste am- 
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OUTLET 
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Figure 3—Venturi scrubber and separator. 
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DESIGN AND PLANT PERFORMANCE 


Dr: SIGN Dat A 
Forecast Load. 


Maximum Load..... 


PERFORMANCE DATA 
(NO expressed as NO2) 


Nov. 1956.. 

Dec. 1956... 

Jan. 1957... 
Gas Rate, cu. 
NO, into absorber, lb. 
NO into absorber, lb. 


NO, in absorber off-gas, Ib. 
NO in absorber off-gas, Ib. 
NO, in Venturi off-gas, Ib. 
NO in Venturi off-gas, Ib. 


Final concentration of outlet gas, p.p.m. 


oxides of nitrogen 


Colume Dia im., Bis ess 
Flow Area, sq. ft.. 
Size of Rings, in.. 


Voidage 


Area of Packing, — sq.ft./cu.ft. 
Gas Rate, cu.ft./min.. . 
Gas Velocity 

through Packing, ft./sec 
Gas Packing Factor.......... 
Liquid Rate, U.S.gal./min. 
H.T.U., ft.. 
* 


ae P: 


ft. 


hr.. 


hr.. 


hr... 


hr.. 


OE: < 


‘hr... 


MPs 2 


COMPARISON OF 


Cases 1, 2, and 3 refer to Aug., 


TABLE 9 


icking physical properties Saini hed by supplier. 


ABSORPTION EFFICIENCIES 


n, lb./hr Out, Ib./hr. Absorption Eff. ‘ 
NO» NO NO, NO and NO» 
78 8 28 65 
480 57 137 65 
91 24 74 
47 11 75 
76 12 84 
TABLE 10 
PLANT PERFORMANCE DATA 
Aug. 1957 _ Sept. 1957 Oct. 1957 
340,000 340, 000 300 ,000 
105 60 75 
| 
25 14 18.4 
18.5 | a7 4 12.5 
1.7 1.6 | 4.4 
2.5 20 | 22 
0:3: | 0.7 0.5 
105 102 115 
TABLE 11 
PILoT PLANT, DESIGN, AND PERFORMANCE 
Pilot i 
Plant Design Calculated 
Case 1 & 2* Case 3* 
1 8 8 8 
0.785 50. 50.2 50. 
1 2 2 2 
0.80 ** 0.74** 0.74** 0.74°° 
46.7 2o.7 ** 4 = pu 
60 6,850 5,670 5,000 
1.6 3.1 Z.00 2,29 
2 7 *** 3.25*** 3.25*** 3.25*** 
a2 250 250 250 
20 34.3 32.2 31.4 
and Oct. 1957 respectively. 


cg 


monia-caustic solution, the scrubbing 
medium has been used on a continu- 
ous bleed-off and make-up basis. The 
average exhaustion of the liquor to 
the sewer is 63%, with exhaustion 
values ranging from 40 to 100%, the 
actual value depending on the avail- 
ability of the waste liquor which was 
completely passed through the ab- 
sorber. Vaporized ammonia is added 
automatically to control the pH. 

More recent data have become 
available. These have resulted from 
the employment of better analytical 
methods and are recorded in Table 
10. 


It is of some interest to compare 
the calculated values of H.T.U. with 
those actually obtained, based on 
total NO and NO, absorption. Table 
11 provides this comparison. The 
“Design” and “Experimental-Calcu- 
lated” values were computed accord- 
ing to reference (23), for a gas-film 
rogryeee absorption process (5 and 

. A satisfactory explanation cannot 
be offered for the lower H.T.U. 
values and improved performances 
obtained by the experimental tests. 

The 
to be 
heavy 


Venturi scrubber has proved 
satisfactory in reducing the 
white ag of ammonium 
nitrate and nitrite. A demonstration 
of its effectiveness can readily be 
observed by noting the change in the 
appearance of the stack gases after 


. VALUES 
Experimental 


Observed 


Case 1* Case 2* Case 3* 


***Gas packing factors obtained from reference (23) and corrected by an experience factor of 20% for the 2 in. packing. 


10 


The Canadian Journal of Chemical Engineering, February, 1958 





wate! 
redi ] 
near: 
chars 
densi 
remo 
nitrit 
NO, 
remo 
only 


Haza 


In 
by a 
and 
nitrit 
is ar 
detor 
temp 
amm 
slow] 
hood 
liquo 
To 
every 
eratic 
depo: 
trite. 
and 
wher 
the < 

Th 
inspe 
caust 
over 


The 





‘ 
water is added to the scrubber. The 
reddish-brown color of the plume 
nearly disappears and the white dis- 
charge is reduced in volume and 
density. According to Table 10, the 
removal of ammonium nitrate and 
nitrite solids, and gaseous NO and 
NO,, ranges from 50 to 85%. The 
removal of solids, based on nitrate 
only ranges from 92 to 97%. 


Hazard 

In the removal of NO and NO, 
by ammonia liquor, the deposition 
and detonation of solid ammonium 
nitrite are a potential danger as this 
is an unstable compound and will 
detonate especially with increase of 
temperature. Aqueous solutions of 
ammonium nitrite decompose very 
slowly and there should be no likeli- 
hood of decomposition of the weak 
liquor in the fume absorption system. 
To maintain a high level of safety, 
every precaution is taken during op- 
eration to prevent the formation and 
deposition of solid ammonium _ni- 
trite. The surface of the scrubbing 
and exhaust systems are kept wet, 
when fume is not passing through 
the absorption unit. 

The equipment was _ thoroughly 
inspected after the waste ammonia- 
caustic liquor had been added for 
over a year. No visible deposits of 


ammonium nitrite or nitrate were 
found in the equipment including 
absorber, Venturi scrubber, separator 
and associated piping. The removal 
of the top of the absorber disclosed 
a yellowish-green colored droplet 
formation which on analysis showed 
no ammonia although 4% nitrite and 
0.5% nitrate. W ashings from the 
base of the Venturi scrubber and the 
separator showed traces of nitrite, 
5-10 p.p.m. 


Inspection 
periodic basis. 


will continue on a 


Conclusions 


In keeping with the policy of the 
Company, considerable time and ex- 
pense were directed towards the 
problem of eliminating a _ potential 
nuisance, the release to atmosphere 
of oxides of nitrogen which are cor- 
rosive in nature and unsightly with 
their intense reddish-brown color. A 
satisfactory method has been de- 
veloped through a literature search, 
a pilot plant ‘study and the design 
and operation of a fume absorption 
system whose efficiency has been 
improved by plant investigations. By 
adopting the methods presented in 
this paper it is now possible to re- 
move the oxides of nitrogen from 
the effluent gases leaving a nitric acid 
plant. 
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Heat Transfer from Column Wall 
to Bed in Spouted, Fluidized 
and Packed Systems’ 


J. KLASSEN? and P. E. GISHLER® 


Heat transfer data have been obtained for heat flowing from a 
steam jacketed column to a packed, fluidized and spouted bed. The data 
for packed and fluidized beds are compared with generalized correla- 
tions previously proposed. In the case of spouted beds, the heat transfer 
coefficient increases sharply at the onset of spouting and passes through 
a maximum as the gas velocity is increased. The maximum heat transfer 


coefficient is from 25-75% lower for a spouted bed than for a fluidized 
bed under similar flow conditions. Rice, wheat and Ottawa sand were 


used as bed materials. 


oo spouted bed method of fluid 
solid contacting has recently been 
described in detail by Mathur and 

Gishler “). Briefly, a ‘spouted bed is 
obtained by introducing the fluid as a 
jet at the bottom of a coarse granular 
bed so as to produce a core of dis- 
persed particles which rises up the 
centre of the bed and culminates in 
a fountain at the top. The particles 
fall out of the fountain or spout and 
form a dense phase annulus which 
slowly travels down towards the 
bottom of the column where the 
cycle is repeated. The annulus be- 
haves like a moving packed bed since 
there is little relative sidew: ay motion 
of the individual particles. 

The spouted bed technique has 
been successfuly applied to the drying 
of wheat ‘2) and shows promise of 
being suitable for other fluid-solid 
contacting operations. Essential pre- 
requisites for an evaluation of this 
technique in connection with any 
particular possible application include 
a knowledge of the flow characteris- 
tics, heat transfer rates and mass 
transfer properties of spouted systems. 
Much work has already been done 


on the first phase“), and it is the 
aim of the present paper to present 
the results of heat transfer measure- 
ments in spouted Sy stems aS a contri- 
bution to the second field of spouted 
bed behaviour. In particular our con- 
cern is with heat flow from the 
column wall to the bed. In order to 
compare the spouted bed with other 
fluid-solid contacting operations, cor- 
responding data were taken for 
packed and fluidized systems. 


Apparatus 


In almost every case the spouted, 
packed and fluidized beds were con- 
tained in a 1 ft. diameter column 8 ft. 
high. A few runs on a 6 in. diameter 
unit are reported; this apparatus will 
be described later. 

Most of the larger column consisted 
of bell-and- ~spigot type sections of 
“Pyrex” piping mounted on a 60 
steel cone. The air inlet at the bottom 
of the cone was 3 in. in diameter, the 
same size as the air line leading to the 
unit. The inlet diameter was varied 
by inserting screened orifices of dif- 
ferent sizes between the two halves of 
the connecting flange. A straightening 
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section of 4 in. copper tubes 10 in. 
long preceded the air inlet to the 
column. Air flow rates above 60 CFM 
were measured with an orifice meter 
while a rotameter was used for the 
lower flows. 

Heat was transferred from a steam 
jacket through the wall of the column 
into the bed. The heat transfer sur- 
face consisted of an 18 in. long steel 
section, held between two of the 
“Pyrex” sections, 18 in. above the top 
of the cone or 26 in. above the air 
inlet. The adjacent “Pyrex” sections 
minimized longitudinal conduction 
and ensured that the heat reached the 
bed via the steel wall. The heat trans- 
fer section is depicted in Figure 1. It 










PYREX SEOTION —PACKING 


STEAM STEAM JACKET 


THERMOCOUPLES (¢ 
TRAVE LLING 

THE RMOCOUPLE 
PACKING 


Ss TO STEAM 
TRAP 


PYREX SECTION 
Figure 1—Heat transfer section. 


consisted of a 12 in. diameter cylinder 
of 3 in. thick steel surrounded by a 
133 in. diameter steam jacket. The 
jacket was terminated 6 in. from the 
bottom of the section in order to 
leave room for the bell- and- “spigot 
connection with the “Pyrex” section. 
A steel bell at the top of the section 
facilitated connection with the upper 
“Pyrex” portion of the column. The 
outside of this heat transfer section 
was insulated with 43 in. of asbestos 
sheet insulation. 
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Che 6 in. diameter unit was similar TABLE 1 
» the 12 in. unit except that it was PHYSICAL PROPERTIES OF THE VARIOUS MATERIALS 


only 33 ft. high with the heat transfer ——- -— = ——SS = SSLEFEFFFESSSEE 


section 6 in. above the air inlet. Bulk Aiea Particle 
Flanged connections were used with Material Densities | Densities Sizes 
the steam jacket extending over the Ib./cu.Et. | Ib./cu.ft. in, 
whole 6 in. length of the heat transfer Celene ena (20 ny 98.1 165 0.0232 
section. 

Temperature measurements were Alundum spheres... .. ars 125 216 0.137 

ade with copper-constantan thermo- 7 : . ‘s 
ne PI ile ns. 56. 92.5 0.187 x 0.093 x 0.066 
couples. Two were embedded 6 in. 
apart, one below the other, in the wall Wheat... 46.4 86.6 0.206 x 0.117 x 0.105 
of the heat transfer section to read 
the wall temperature. A third thermo- 
couple was installed on a level be- 
tween the first two in such a way that termined with a micrometer. An through the unit and the steam flow 


it could be moved horizontally across 
a column diameter to read radial 
temperature profiles. The thermo- 
couple leads passed out through the 
steam jacket and then to a Rubicon 
potentiometer and ice-water system 
where the thermal emf’s' were 
measured. 

The heat transfer rate was calcu- 
lated from the enthalpy of the in- 
coming and outgoing steam and the 
steam flow rate. The initial enthalpy 
was determined from the temperature 
and pressure of the incoming steam, 
which had been superheated 5 to 
10°F. with a 100 watt electric heater 
inserted in the steam line. After leav- 
ing the steam jacket and passing 
through a float-type steam trap, the 
condensate was led to a calorimeter 
for determination of the final en- 
thalpy and flow rate. 


Experimental and Calculations 


A. Physical Properties of the 
Solids 


Before any heat transfer runs were 
conducted some of the physical pro- 
perties of the various materials used 
were measured. These properties in- 
cluded bulk densities, absolute par- 
ticle densities and particle size. Table 
| lists these properties for each of the 
materials. 

The bulk densities were determined 
by weighing the quantity of material 
needed to fill a 2 litre cylinder. The 
absolute densities were determined by 
several methods. For the Ottawa sand 
a water displacement pycnometer was 
used. The values for rice and wheat 
were measured with an absolute alco- 
hol displacement pycnometer. In the 
case of the alundum spheres 25 par- 
ticles were each weighed and the 
diameter measured with a micrometer. 

The particle size of the Ottawa 
sand (20-30 mesh) was taken as the 
average of the openings of the two 
screen sizes. In the case of the wheat 
and rice two of the dimensions of 
each of 25 particles were measured 
with a travelling microscope. The 
third and smallest dimension was de- 


average based on the 25 particles is 
reported for each dimension. The 
diameters of the spherical alundum 
particles were read from a micro- 
meter. Again 25 particles were taken 
as the sample. 


B. Heat Transfer Coefficient 
The heat transfer coefficient was 
calculated from the observations by 
means of the equation: 
q 


h = .(1) 
: AAt ( 


Heat transfer rate and temperature 
difference were obtained as follows: 


C. Heat Transfer Rate 

The heat transfer rate, q, was 
taken as the difference between the 
heat given up by the condensing 
steam ‘and the heat losses through the 
insulation and steam trap. The former 
quantity was based on the enthalpy 
change of the steam in _ passing 
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rate. The pressure and temperature of 
the incoming steam, which was super- 
heated 5 to 10°F., together with the 
data of the steam tables, permitted 
a determination of the initial steam 
enthalpy. 

A calorimeter sitting on a balance 
pan and equipped with a Beckmann 
differential thermometer was used to 
measure the enthalpy of the steam 
condensate as well as the steam flow 
rate. The weight of water in the 
calorimeter was about 40 lb., while 
the amount of condensate collected 
per determination was about 2 lb. 
Weighings were accurate to within 
0.1 oz. The calorimeter determination 
was timed with a stop-watch; the time 
varied from a few minutes to over an 
hour. A weight calibration of the 
calorimeter showed that the buoyant 
effects on the thermometer, agitator, 
etc. were negligible while a thermal 
calibration showed that the heat capa- 
city of the unit, based on the weight 
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Figure 2—Vertical temperature profiles. 
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Figure 3—Caleulated and observed values of h for a packed bed of rice. 


TABLE 2 


HEAT TRANSFER IN PACKED BEDs* 


Inlet \ir Flow Rate ts tb At q h 
Diam. : ; ; a 
in. G,lb./(hr. ) F. F. i BTU/hr.) BTU/(hr.) 
(sq. ft.) (sq.ft (°F. ) 
perforated 
grid air 
inlet 239 252 74 156 1530 232 
371 252 75 151 2210 3.3 
196 252 77 146 2870 4.4 
con al air 
inlet 1 456 251 79 147 1820 2.8 
2 182 254 80 159 960 1.4 
2 371 253 80 150 1870 2.8 
2 480 255 81 148 2210 3.4 
2 554 253 81 143 2800 4.4 
3 322 254 83 152 1320 1.9 
3 505 253 76 150 2470 3.9 
* Material: rice, a5 2 ft, 


bed depth: 4 ft., fluid: air. 





of water, varied with the time of vertical surfaces, recommended by 
operation. This variation was des- McAdams ‘*) was used. 
cribed by 


: heL Lp;*g3,A Cour Vy 
heat capacity aes 0.13 ( Pp aI ‘)( ph )| 4 
1 + 0.00782 (time in minutes) °-75 ks bs* ky 


The heat losses, i.e. the heat given 
up by the steam which did not reach 
the bed via the heat transfer surface, 


For the conditions existing in the 12 
in. diameter column 


i.e. oe 
were determined from runs with free e. fluid - 

° ° € . 
convection in the empty column. A At = 89°F. 
correlation for free convection at A = 4.48 sq. ft. 


the above equation gives h, = 0.87 
BTU/hr. sq. ft.°F. 
Heat lost by convection 
q. = h,A At 
= 350 BTU/hr. 

Heat is also lost from the heat 
transfer section by radiation. The co- 
efficient of heat transfer by radiation 
is given by the following equation: 


h, = 


t.+460\' tp+460\' 
oataxexl (EM) 4 (SHY) 
100 100 


ts— tb 


It is assumed that the heat lost by 
radiation is equivalent to that lost by 
two discs equal in area to the cross- 
section of the column and having the 
same temperature (252°F.) and emis- 
sivity (0.66) as the heat transfer sur- 
face. Heat from the top disc is then 
lost to that part of the column which 
is at the bed or warm air temperature. 
That from the lower disc is lost to 
the column which is at ambient tem- 
perature. Substituting these conditions 
1.e. bed temperature 150° and ambient 
temperature 80°F., the coefficients of 
heat transfer are 1.32 and 1.14 
BTU/hr.sq.ft.°F. for the top and 
lower ends respectively. 

Cross-sectional area of column = 
0.785 sq. ft. 

Heat lost from top end, 

q = h,AAt 

= 1.32 x 0.785 x (252-150) 

= 106 BTU/hr. 

Heat lost from lower end 

= 1.14 x 0.785 x (252-80) 

= 154 BTU/hr. 

Hence total heat passing through heat 
transfer section and into the column, 

a = 350 + 106 + 154 

= 610 BTU/hr. 
Total heat given up by steam = 1440 
BTU/hr. Heat losses in steam pipes, 
jacket, insulation, etc. are, 
1440-610 = 830 BTU/hr. 

Variation of steam jacket temper- 
ature had little effect on the heat 
transfer coefficient. Most runs were 
therefore carried out at the same 
temperature and the heat loss for all 
runs was constant. A correction of 
830 BTU/hr. has been applied in each 
case. 

D. Temperature Difference 

The temperature difference, At, 
was calculated from the surface tem- 
peratures measured by the two 
thermocouples embedded in the 
column wall and the average bed tem- 
perature as measured by the travel- 
ling thermocouple. Measurements of 
the vertical temperature profile were 
made to see if a bed temperature 
reading at one level (i.e., that of the 
radially -travelling thermocouple ) 
would be sufficient for calculating At. 
This profile together with the cor- 
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responding wall temperatures to be 
discussed subsequently, is given in 
Figure 2. The uniformity of the bed 
temperature justifies the use of the 
bed temperatures read on the one 
level. 

Fhe wall temperature profile is 
complicated by the fact that the 
steam jacket extends only a part of 
the way down the heat transfer sec- 
tion. As would be expected and as 
was born out by the readings of the 
two embedded thermocouples the 
wall temperature is constant over that 
portion of the surface backed by the 
steam jacket. A temperature gradient 
is to be expected, however, over the 
rest of the surface. This gradient a 
can be formulated mathematically by bicep : ee _| 
solving the differential equation for 9 200 400 600 800 1000 ~—Ss«1200 1400. 1600 1800 
steady state conduction in rods and G, MASS VELOCITY, LB./(HR)(SQFT.) 


other projections with the aid of ap- 
propriate boundary conditions. This Figure 4—h vs G curves for fluidized beds of rice, wheat, sand, and alundum 
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has been done in the appendix. The spheres. 
resulting wall temperature profile has 
been included in Figure 2 together 
with the true At which, according to 
the appendix, can be ca!culated as 


) 
As + : tanh ml a 
a m i. — &). (Se) PABLE 3 
A HEAT TRANSFER IN FLutIpIzED BEDs* 
= @(ts — th)... » «(ODD 
The introduction of equation (5b) - - Air Flow Rate ‘ tb At q h 
into equation (1) gives Materia : . ; dee os 
(1) g Glb./(hr.) | °F. | °F. | °F. | BTU/hr. | BTU/(hr.) 
= (h@) A (ts tn) <3) (sq. ft.) (sq.ft )(°F.) 
Equations (5a) and (5b) show that 4 Ottawa sand 126 252 77 159 970 1.4 ‘ 


is a function of h. By choosing ag 


e »? Ss 75 S52 g 
of h and calculating values of ¢ : 408 Zo! = 153 1800 2.6 


graph of @ vs h or, as is more con- 372 250 185 46.7 8220 39.3 
venient, a graph of hd vs @ can be 
plotted. For each run a value of h¢ | 486 249 178 51.6 9010 39.2 
can be calculated from equation (3) | 991 249 134 84.7 10050 5 6 
The graph referred to above gives a 
corresponding value of ¢ which al- Alundum spheres | 628 251 | 77 | 143 3300 $4 
lows a calculation of the heat transfer | 1724 48 +43 we 14040 a 
coefficient to be made h = (hd@)/¢. si ve . 7 ; = 
Rice 239 252 74 156 1530 aon 
Results 
371 252 75 5 2210 in 
A. Packed Beds " E " = . ' 
The majority of the packed bed 496 | 252 | 77 145 2870 4.4 
heat transfer data were taken for beds se nee 
557 250 150 74.3 


fat tran: 6960 20.9 
of rice. Table 2 and Figure 3 show iid 


the results. The fact that the data 609 251 143 89.2 | 2740 4.7 
obtained with a perforated grid coin- “3 . 
cide with the data obtained with a | ae | = 136 | of. oH 42.9 
spouting-type air inlet indicates that 736 | 251 135 85.7 0670 25.2 
the air distribution must have been 
similar in both cases. 1060 | 250 | 123 94.3 10310 | 22.3 

:xtrapolations of correlations of Chines | 492 | 957 | 74 | 447 2960 4.5 
several workers in the field of wall to 
packed bed heat transfer have been 683 252 75 144 3440 5.2 
included in Figure 3 (45.5.7), The wide 819 | 050 128 90.5 8220 10.3 
differences are most striking. This , vr - “i . oa 
seems to be due to the extrapolation 991 250 i | | |69S:% 8750 20.4 
of the correlations, which were de- 

1219 250 117 99.2 9450 ania 


veloped from data in columns under 
6 in. in diameter, to a column = 
diameter of 12 in. There is an indica- * De: 1ft., bed depth: 4 ft. 
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Figure 6—h vs G curves for spouted beds of sand, rice and wheat. 


tion that generally the column diam- 
eter has only a minor effect on the 
heat transfer because if the correla- 
tions are modified by ignoring the 
factors involving the diameter they 
tend to coincide more closely. 


B. Fluidized Beds 


Fluidized bed heat transfer data 
were obtained with sand, alundum 
spheres, rice, and wheat as_ bed 
materials. These materials, particu- 
larly the latter three, fluidized poorly 
because of severe bumping. The data 
are contained in Table 3 and plotted 


16 


in Figure 4. The general shape of the 
curves with a maximum value of the 
heat transfer coefficient is in agree- 
ment with previous work (8) as is the 
trend of heat transfer with partic'e 
size, Le., the larger the particle the 
lower the maximum heat transfer co- 
efficient. The maximum is the resul- 
tant effect of two opposing mechan- 
isms. With increasing air flow rate 
the turbulence in the bed increases 
tending to increase h. The increasing 
flow rate also causes an expansion of 
the bed with a corresponding decrease 
in interstitial air velocity. At low air 


flows the first effect predominates 
while at high flows the latter effect 
governs the behaviour. Somewhere in 
the intermediate region the two op- 
posing effects must produce a maxi- 
mum. 

The experimental curve for rice is 
shown in Figure 5 along with curves 
computed from various proposed 
formulas (9:10,11,12,13,14,15), Tt is evident 
that only a few of these formulas came 
close to predicting the heat transfer 
coefficients actually obtained in the 
twelve inch unit. Only correlating 
formulas by Van Heerden et al (!!) 
and Toomey, R. D. I. and Johnstone, 
H. F. “5) take into account the fact 
that the heat transfer rises suddenly 
at incipient fluidization, but even 
these do not predict the shape of the 
heat transfer curve over a velocity 
range. Clearly, a generalizee heat 
transfer correlation is not yet avail- 
able (as already noted by Van 
Heerden ‘'®)) for predicting fluidized 
bed heat transfer coefficients over a 
gas velocity range. 


C. Air Spouted Beds 


The data on heat transfer in 
spouted beds are given in Tables 4 
and 5. Some of the results of Table 4 
have been plotted in Figure 6. The 
general shape of the heat transfer co- 
efficient vs. mass velocity curve is very 
similar to the curve for a fluidized 
bed. In both cases there is a sharp 
increase at the point of transition 
from a packed to a moving bed. The 
heat transfer coefficient rises to a 
maximum after which it falls off with 
increasing velocity. The coefficient 
for a spouted bed is invariably lower 
than for the corresponding fluidized 
bed. This is to be expected as there 
is greater turbulence of the particles 
in the latter compared with the 
smooth, orderly, laminar-like motion 
of the particles at the wall in the 
former. 


Apparently the smaller the particles 
the lower the heat transfer coefficient 
for the spouted system and _ the 
greater the difference in heat transfer 
rates between the spouted and flui- 
dized state. This is evident when the 
data for sand (particle diameter of 
about 0.02 inches) is compared with 
that for wheat and rice (particle 
diameter of about 0.15 in.). The 
larger bed turnover rates associated 
with the spouting of larger particles 
may give rise to the larger heat 
transfer coefficients. This shows that 
the effect of particle diameter on heat 
transfer coefficient in a spouted bed 
is the reverse of the trend for a flui- 
dized bed. Too few materials were 
available to make a more reliable study 
of this variation. 
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TABLE 4 


HEAT TRANSFER IN SPOUTED BEDS 

















Nl l | | | 
Column | Inlet | Bed | Air Flow Rate | ts to Mm | q h 
Material Diam. | Diam. | Height | 

in. in. | in. | G,lb./(hr.)(sq.ft.)| °F. | °F. °F: BTU/hr. | BTU/(hr.)(sq.ft.)(°F.) 

Ottawa sand | 12 0.5 | 4 229 | 251 101 118 | 4920 9.4 

| 12 0.75 4 272 | 250 97 119 | 5250 9.9 

a | 68 | C4 352 251 104 115 | 5040 9.8 

Rice 12 1 | 4 456 251 79 148 1820 2.8 

12 1 | 4 556 248 86 124 7530 13.6 

12 1 | 4 616 | 249 | 93 120 6840 12.8 

12 1 | £ 652 | 248 90 121 6480 11.8 

42 1 |} 4 | 681 250 96 112 6910 13.1 

12 2 | 4 182 | 254 80 159 960 1.3 

12 2 | 4 371 } 254 80 150 1870 2.8 

12 2 | 4 480 | 255 81 148 2210 3.4 

12 2 | 4 554 | 25a- | 81 143 2800 4.3 

12 2 | 4 572 251 | 88 125 7030 12.5 

12 2 | 4 628 253 | 96 17 || 9760 17.1 

12 2 | 4 677 251 | 91 120 | 8570 15.9 

12 2 | 4 720 252 | 90 121 10250 19.2 

ym | 2 | 4 | 794 | 252 91 121 8750 16.1 

12 2 | 4 | 862 | 251 90 121 9750 | 18.1 

12 3 4 | 322 | 254 83 152 1310 | 3.5 

12 3 4 505 | 253 76 149 | 2590 | 3.9 

12 3 4 572 251 97 119 7010 13.2 

12 3 4 733 251 91 7 | 8040 | 14.7 

12 3 4 | 886 251 91 20 | 8600 15.9 

6 1 e | 543 248 80 167 140 1.1 

6 1 / a 606 249 136 112 1510 17.2 

6 1 2 | 650 248 127 120 1890 20.0 

Wheat | 12 1 a | 681 247 100 113 6460 13.2 

| 12 1 4 | 804 248 96 116 | 7030 | 13.5 

12 1 | 4 | 348 248 93 118 | 8260 | 13.7 

| 12 1 | 4 939 248 93 117 7900 15.0 

12 3 | 4 | 322 252 76 150 2190 3.3 

| 12 3 & | 677 249 96 | 131 7490 14.2 

| 12 3 4 | 751 251 | 81 | 140 | 3040 4.9 

12 3 4 | 819 249 | 91 | 119 | 8360 15.6 

The maximum heat transfer coeffi- TABLE 5 


cients at different inlet diameters are 
given in Table 5, for the various 
materials of this study. It can be 
concluded that the inlet diameter had 
no effect, or that the particle break- 
up, which was much more pro- 
nounced for the smaller inlet sizes, 
changed the effective partic'e diam- 
eter so as to counteract any effect of 
the inlet diameter alone. The particle 
break-up was a serious problem in 
this work. Each run required a period 
of one or more hours for the attain- 
ment of steady-state conditions. Even 
a small rate of particle break-up 
could have produced major changes 
in the system after such long periods. 
The materials included in this project 
were chosen after numerous solids 
had been examined. Many were 
quickly ruled out after their exces- 
sively high attrition rates were ob- 
served. The materials which were 
studied were therefore chosen pri- 
mari'y because of their attrition re- 
sistance. 

It should be noted that column 
diameter had no apparent effect on 
the observed heat transfer coefficient. 
Pertinent data obtained with rice 
using the 6 in. and 12 in. column are 
listed in Table 4. 


MaAxiImMuM HEAT TRANSFER COEFFICIENT IN A SPOUTED BED USING 


Material 


Ottawa sand.. 





12 in., bed height: 4 ft. 


Summary of Results 


1. In fluidized beds, heat transfer 
coefficients increase suddenly on the 
onset of fluidization, and pass through 
a maximum value as mass velocity is 
increased, confirming previous data 
obtained in a smaller unit. 


2. In spouted beds, heat transfer 
coefficients follow a similar relation- 
ship with increasing mass velocity, 


The Canadian Journal of Chemical Engineering, February, 1958 


Various Air INLET DIAMETERS* 


Maximum heat transfer 





Inlet Diameter coefficient 
in. BTU/(hr. )(sq.ft. )(°F.) 
\ 9.4 
34 9.9 
1 13.6 
2 19.2 
3 15.9 
1 15.0 
a 15.6 


but the maximum values are smaller 
—about 25 to 75% lower than for a 
fluidized bed under similar conditions. 


3. The maximum heat transfer co- 
efficient in a spouted bed is indepen- 
dent of inlet orifice diameter, and of 
column diameter. There are indica- 
tions that h,,,, increases with bed 
particle diameter rather than decreases 
as normally observed with fluidized 
beds. 
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Appendix 


A. Calculation of Wall Temper- 
ature Profile and an Effective 
At for Heat Transfer 


The wall of the heat transfer sec- 
tion backed by the steam jacket (see 
Figure 1) will be at the condensing 
steam temperature. The portion of 
the wall below the level of the steam 
jacket will contain temperature grad- 
ients. A heat balance around an ele- 
ment of this part of the wall results 
in the differential equation: 

Se nM Wao sccn cal 
dx? 
The section on nomenclature defines 
all the letters and symbols. 

Equation (4) can be solved with 

the boundary conditions: 


T = T,atx =O 
(dT/dx) = Oatx =1 
to give 
T 


= cosh mx — tanh ml sinh mx 
. 


The effective temperature differ- 
ence can be calculated from 


1 


q = hAAt = hA,T; + Ph | T dx 
or 
P 
As + tanh ml 
OS a ad (ts — tp). .(5a) 
A 
Gs hi os 6 sain RGAE (Sb) 


For the 12 in. diameter column 


4.48 sq. ft., a = 0.0315 sq. ft., 
0.5 ft. 


k = 26 BTU/hr. ft. °F. for steel, ref (17) 


~ > 
Nit 


18 


hP hxamx 11.5 
m? = = - = 3.67h ft.-* 
ka 12x 26x 0.0315 
m = 1.917 h!/? ft.-! 
ml = 0.958 h!/2 
P wx 11.5 


Ve — 157 h-1/2 ft 2 
m 12x 1.917 xh! — os 


Substituting in equations (Sa) and 


(5b) 


2.99 + 1.57h~!/? tanh 0.958 h!/? 


4.48 

Nomenclature 

A = area through which heat flows, sq. 
ft. 

A, = heat transfer area backed by steam 
jacket, sq. ft. 

a = cross sectional area through which 
heat is being conducted = mDc x 


(thickness of metal), sq. ft. 


C, = heat capacity, BTU/(Ib.) (°F.). 

D. = column diameter, ft. 

G = mass velocity, lb./(hr.) (sq. ft.) 

g = gravitational acceleration = 4.2 x 
108, ft./(hr.) (hr. ) 

h = heat transfer coefficient for heat 
flow from the column wall to the 
bed, BTU/(hr.) (sq. ft.) (°F.) 

h. = heat transfer coefficient for free 
convection, BTU/(hr.) (sq. ft.) 
(°F.) 


hmax = Maximum heat transfer coefficient, 


BTU/(hr.) (sq. ft.) (°F.) 


h, = heat transfer coefficient for radia- 
tion BTU/(hr.) (sq. ft.) (°F.) 

k = thermal conductivity of a solid, 
BrU/(hr.) Gt.) CF.) 

kp = thermal conductivity of a fluid, 
BTU/(hr.) (ft.) (°F.) 

i = height of heat transfer section, ft. 

l = vertical distance between the bot- 


tom of the steam jacket and the 
bottom of the heat transfer sec- 


tion, ft. 

m «= h P/ka, ft.-? 

P  =circumferential perimeter of heat 
transfer section = 7D,, ft. 

q = heat flow rate from column wall to 


the bed, BTU/hr. 


T = tw — tn, °F. 
Tr, =t. — to, F. 
t = temperature difference between 


column wall and bed, °F. 
tp = bed temperature, °F. 


ts = temperature of condensing steam, 

tw = wall temperature, °F. 

x = vertical distance below the bottom 
of the steam jacket, ft 

Bi = fluid coefficient of thermal expan- 
sion, (°F.)~! 

bf == fluid viscosity, Ib./(ft.) (hr.) 

w = 3.14159 

pi = fluid density, Ib./cu. ft. 

@ = (A, + P/m tanh ml/A) 

€ = emissivity of heat transfer section. 
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Relative Volatility and Enthalpy 
Data for the Systems 
C, Hydrocarbons-Acetone- Water 
Developed from Vapor-Liquid Equilibria’ 


J. E. EWANCHYNA? and C. AMBRIDGE® 


Extractive distillation is widely used in separating C, hydrocarbons. 
A typical extractive solvent for such a separation is an acetone-water 
mixture. Correct design and selection of optimum operating conditions 
for extractive distillation towers are ensured by a knowledge of the vapor- 


liquid equilibria involved. 


Pseudo-binary systems of each of the C, hydrocarbons, n-butane, 
isobutane, butene-l, and cis-butene-2, and acetone-water solvent were 
investigated at temperatures of 100°, 150°, and 200°F., at water concen- 
trations of 5, 11, and 17 wt. %, (C, hydrocarbon-free) and at hydrocarbon 
concentrations up to 28 mole %. The results were interpreted as activity 
coefficients, and showed good consistency. The binary Margules three- 
suffix and van Laar two-suffix equations correlate the measured activity 
coefficients adequately as a function of concentration. 

Dependence of activity coefficients on temperature was used to 
evaluate the differential heats of solution of the hydrocarbons in the 
acetone-water solvent. The results were also used to derive partial liquid 


enthalpy-concentration curves. 


Ternary Margules three-suffix and van Laar two-suffix equations 
gave a satisfactory correlation of activity coefficients for pseudo-ternary 
systems. The ternary Margules equations were used to calculate relative 
volatilities, which are presented graphically as functions of temperature 


and composition. 


i ensure correct design and eco- 
nomic operation of extractive 
distillation towers used in separations 
of C, hydrocarbons, knowledge is 
required of the vapor-liquid equilibria 
that are involved. Such information is 
needed particularly for multicom- 
ponent systems which deviate con- 
siderably from ideality, because the 
relative volatilities of the hydro- 
carbons in the extractive solvent de- 
pend on the total concentration of 
hydrocarbon in the solvent and the 
relative concentrations with respect 
to each other. Two widely used ex- 
tracting solvents are acetone-water 
mixtures and furfural-water mixtures. 

Accurate determinations have been 
reported of vapor-liquid equilibrium 
data on C, hydrocarbons with dry 
furfural (1: 2) and with furfural-water 


mixtures (34), These binary, ternary 
and quaternary mixtures of non-ideal 
components were successfully treated 
with the Margules three- suffix and 
van Laar two-suffix equations which 
were derived by Carlson and Colburn 
(5) and Wohl (6). It was shown (3,4) 
that ternary mixtures of hydrocarbon 
with furfural and water, and quater- 
nary mixtures of two hydrocarbons 
with furfural and water, could be 
simplified to pseudo- -binary and 
pseudo-ternary mixtures respectively. 
This could be done by treating the 
solvent, furfural-water, as one pseudo- 
component as long as the water con- 
tent of the solvent remained constant. 
Such treatment of data greatly simpli- 
fied expressions for calculating rela- 
tive volatility, reduced the number 
of systems that needed to be studied, 


1Manuscript received November 18, 1957. 
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and enabled easier graphical represen- 
tation of data. 

The object of this investigation was 
to measure vapor-liquid equilibrium 
data for C, hydrocarbon-acetone- 
water systems for derivation of rela- 
tive volatility and enthalpy data. 
Pseudo- binary systems of each of the 
hydrocarbons, ‘n-butane, isobutane, 
butene-1, and cis-butene-2, and the 
solvent, acetone-water, were investi- 
gated at temperatures of 100°, 150°, 
and 200°F., at water concentrations 
of 5, 11, and 17 wt. %, (C, hydro- 
carbon-free) and at hydrocarbon 
concentrations up to 28 mole %. 
The studies were confined to a range 
of approximately 0 to 28 mole % 
hydrocarbon to avoid phase separa- 
tion. For all calculations the acetone- 
water solvent was treated as a single 
pseudo-component. A pseudo-ternary 
system of n-butane-butene-1-acetone- 
water was investigated at 100° and 
150°F., with a water concentration of 
11 wt. % (C, hydrocarbon-free). 


Materials 

Pure grade acetone containing up 
to 0.4 wt. % water was used. The 
solvent mixtures were made up by 
adding the calculated weight of water, 
followed by readjustment w ith water 
or acetone to give the required water 
content as determined by analysis. 

The hydrocarbons were obtained 
from Phillips Petroleum Company 
and were the pure grade, specified 
over 99% pure. 


Description of apparatus 

A flow-type of apparatus was used 
in which a vapor stream flowed 
through a liquid mixture at constant 
temperature and pressure until equi- 
librium was attained. A schematic 
flow diagram is shown in Figure 1. 
This apparatus is a slightly modified 
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Figure 1—Equilibrium apparatus. 
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Figure 2—Stripping apparatus. 


version of that used at the University 
of Delaware “), 

Liquid hydrocarbon stored in the 
supply tank was heated to keep its 
vapor pressure at a level sufficient to 
maintain a flow through the apparatus. 
Liquid hydrocarbon was metered 
through a rotameter and vaporized in 
the flash vaporizer. The hydrocarbon 
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vapor passed into the first saturator 
chamber through a heater, which 
maintained the temperature of the 
vapor a few degrees above that of the 
saturator. The vapor next passed 
through a second saturator into the 
equilibrium chamber, both vessels 
being maintained at the same temper- 
ature by thermostatic control. The 










second saturator and equilibrium 
chamber, and was normally main- 
tained a few degrees higher. 

Both saturators contained solvent 
to ensure that the vapor stream enter- 
ing the equilibrium chamber was 
saturated with solvent vapor and 
would not change in composition as it 
passed through the solvent in the 
equilibrium chamber. Use of a second 
saturator was the major modification 
of the apparatus from that used at the 
University of Delaware. The com- 
position of the liquid in the first 
saturator alters slightly due to the 
fact that the acetone/water ratio in 
the vapor differs from. that in the 
liquid; without the second saturator 
this would affect the composition of 
the vapor entering the equilibrium 
chamber. In practice the effect was 
negligible at 100° and 150°F., but was 
appreciable at 200°F. With the second 
saturator the composition of the 
vapor entering the equilibrium cham- 
ber was constant at all three temper- 
atures. 

The vapor stream leaving the equi- 
librium chamber was condensed and 
collected in a receiver. Constant pres- 
sure in the equilibrium chamber was 
maintained by carbon dioxide back 
pressure. Diffusion of carbon dioxide 
to the equilibrium chamber was pre- 
vented by the flow of hydrocarbon 
through the condenser. 
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Sampling procedure 

Sample bombs were made of brass 
or aluminum, and were fitted at each 
end with a 4 in. Hoke valve. The 
size of bomb ranged from 2 ml. to 50 
ml., so that samples of any required 
volume could be taken, depending on 
the composition of the sample and 
the analysis required. 

In taking liquid samples, a weighed 
evacuated bomb was securely con- 
nected to the sample line at the 
bottom of the equilibrium chamber. 
The valves were opened to allow 
liquid to fill the bomb, then the bomb 
outlet valve was cracked open and 
the bomb was slowly flushed with at 
least four times its volume of liquid. 
The valves were closed, and the bomb 
was disconnected, dried and weighed 
to determine the weight of sample. 

In taking vapor samples, it was 
necessary to condense a sample of 
vapor into a bomb in order to obtain 
sufficient sample for analysis. It was 
important to prevent condensation 
occurring in the line before the 
sample point, and to ensure that 
complete condensation took place as 
the sample entered the bomb. The 
sample point was located as close as 
possible to the equilibrium chamber, 
and the temperature of the vapor was 
maintained 5°F. above that of the 
chamber by means of a heating tape 
wrapped around the line. A weighed 
evacuated bomb was connected to the 
sample point, the space in the union 
connection was evacuated, and the 
bomb and connection were carefully 
tested for leaks. The lower end of 
the bomb was cooled with liquid 
nitrogen, while the line up to the 


TABLE 1 


EXPERIMENTAL COMPRESSIBILITY FACTORS AT APPROXIMATELY 75°F. AND 
Za 0 ATMOSPHE RE 














Mean Number | Standard | Literature | 
Hydrocarbon Experimental | of Deviation | Value Reference 
Determination | Determinations | | 
n-Butane 0.9732 7 0.0025 0.9728 1 
Butene-1 | 0.9711 | 9 0.0046 | 0.9703 | 1 
Isobutane | | 0.9754 | 1 
cis-Butene-2 | 0.9741 | 4 


valve was maintained at a high tem- 
perature. The valve was cracked 
open, and a sample slowly condensed 
into the bomb. The sample was taken 
slowly in order to maintain a flow of 
vapor through the condenser to pre- 
vent carbon dioxide from leaking 
back to the sample point. 

A run was usually completed in less 
than 8 hours. The level of liquid in 
the equilibrium chamber was ob- 
served to rise as hydrocarbon was 
dissolved, and equilibrium was indi- 
cated when this level remained steady. 
This usually occurred about two 
hours from the start of the run, and 
sampling was then started. Three sets 
of samples were taken for each run. 
Each set consisted of two liquid and 
two vapor samples, because in each 
phase one sample was required for 
hydrocarbon and acetone analysis, 
and one for water analysis. Normally 
two sets of samples were analy zed, 
and if the agreement was good the 
third set was discarded. For some 


TABLE 2 





runs only one set of results has been 
recorded; in these cases one or more 
samples were accidentally spoiled and 
the complete analysis was not avail- 
able, but in no case was a result 
accepted unless the hydrocarbon 
analyses on two successive liquid 
samples showed that equilibrium had 
been attained. 


Analytical procedure 

The scheme of analysis for 
hydrocarbon and acetone was to flush 
the sample from the bomb with 
carbon dioxide, remove acetone by 
passing the stream through an aque- 
ous solution of hydroxylamine hydro- 
chloride, and measure the volume of 
hydrocarbon in a gas buret after re- 
moving the carbon dioxide by absorp- 
tion in a concentrated caustic solution. 
Figure 2 is a schematic diagram of 
the apparatus. Acetone was deter- 
mined by titrating the hydroxylamine 
hy drochloride solution with standard 
alkali solution. Water was determined 
on a separate sample by titration with 


PsEUDO-BINARY VAPOR-LIQUID EQUILIBRIUM DATA FOR n-BUTANE AND BUTENE-1 IN ACETONE-WATER SOLVENT CONTAINING 
Ls 4 W1 T. % W. ATE R (Ce Hy DROCARBON- FREE) 














| Vapor Analysis (Mole %) Liquid Analysis (Mole %) | Activity 

Total se aa hdl ee einige es Correction | Coefficient 
Run No. Pressure Factor of 

Mm. Hg. | Acetone Water Hydrocarbon | Acetone Water Hydrocarbon z Hydrocarbon 
n-Butane, Temperature 150.3°F. 
100 ie tl CORES OY CSC 46.8 82.3 | 4:5 | 3.20 | 1.128 | $.35 
1769 arG | S68 | 47.3 81.9 14.8 3.28 1.128 Sao 
99 2787 | yn ee. | 68.7 76.6 14.4 9.00 1.092 | 4.24 
2784 206.6 | 3.6 | 68.8 76.5 14.6 8.94 1.092 4.27 
116 3534 20.1 So. 76.8 71.4 13.4 15.21 1.065 3.47 
105 4303 16.9 2.4 80.9 62.7 | 11.6 25.03 1.039 2.56 
4303 16.2 2.0 81.8 62.0 2.1 25.92 | 1.039 2.58 
Butene-1, Temperature 150.4°F. 

113 Moe 6) ae) CL Se | 46.4 81.8 14.0 | 4.20 1.151 3.45 
1764 St.2 | 5.5 46.6 81.6 14.2 4.17 1.151 3.49 
114 2788 26.3 | 4.0 69.7 74.6 14.0 11.42 1.117 2.92 
2788 26.1 | Ge 69.6 74.9 13.8 11.34 1.117 2.94 
111 | $535 5 a 75.8 68.7 13.4 17.92 | 1.092 2.51 
| 3535 18.7 3.0 77.8 68.8 12.9 18.34 |} 1.092 2.92 
108 4303 15.2 2:0 82.1 60.4 11.7 27.92 | 1.068 2.07 
15.9 ana 7 | 60.6 Rid 28.09 1.068 2.05 


| 4303 


| 81. 
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Figure 3 — Activity coefficient values for n-butane in acetone-water solvent 
containing 11.0 Wt. % water (C, Hydrocarbon-Free) showing effect of varying 
temperature. Solid curves are predicted by binary Margules equation. Dotted 


curves are predicted by van Laar equation. 
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butene-| in acetone-water solvent 

containing 11.0 Wt. @% water (C, Hydrocarbon-Free) showing effect of varying 

temperature, Solid curves are predicted by binary Margules equation. Dotted 
curves are predicted by van Laar equation. 


Figure 4 — Activity coefficient values for 


Karl Fischer reagent. 


Analysis for hydrocarbon 


The closed sample bomb was 
placed in position, and carbon di- 
oxide was passed through the appa- 
ratus to remove all air. The carbon 
dioxide was absorbed in the 25 wt. % 
solution of potassium hydroxide con- 
tained in the gas collection buret, and 
any impurity was collected in the 
buret. The flow of carbon dioxide was 
continued until successive metered 
volumes of carbon dioxide gave the 
same volume of impurity. The maxi- 
mum amount of impurity was of the 
order of 0.1 ml. per litre of carbon 
dioxide. 

After expelling the collected gas 
from the buret, the bomb was opened 
slowly, and the released hydrocarbon 
vapor passed through the hydroxyla- 
mine hydrochloride solution into the 
buret. A metered volume of carbon 
dioxide was passed through the appar- 
ratus to remove all hydrocarbon and 
acetone from the bomb and to strip 
any dissolved hydrocarbon from the 
hydroxylamine solution. Once more 
the flow of carbon dioxide was con- 
tinued until successive volumes of gas 
gave the same volume of impurity. 
The collected gas was transferred to 
the measuring buret containing fresh 
potassium hydroxide solution, where 
its volume was measured and the 
temperature noted. Solubility of 
hydrocarbon in potassium hy droxide 
solution was negligible ‘!). The vol- 
ume was corrected for the amount of 
impurity in the carbon dioxide; nor- 
mally the value of the correction was 
less than 1 ml. 


The weight of hydrocarbon was 
calculated from Equation 1, 


pVM 
m = pm, vars (1) 

Cm 
The partial pressure of the measured 
volume of hydrocarbon was obtained 
by subtracting the water vapor 
pressure of the potassium hydroxide 
solution “) from the corrected baro- 
metric pressure. Compressibility fac- 
tors for n-butane and butene-l were 
determined experimentally by using 
known weights of pure hy drocarbon 
as samples. Good agreement was ob- 
tained with the values determined by 
previous 


workers ‘'). For isobutane 


and cis-butene-2 the literature values 
'), Compressibility fac- 


tors are given in Table 


were used 


In the case of hy drocarbon mix- 
methods 
determine the ratio of components in 
the mixture. For Runs 20-32 the gas 
sample collected in the buret was 
analyzed by infra-red absorption, and 


tures, two were used to 
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:olecular weights were used in the 
gas law equation. For runs 74-80, 84, 
119 and 124, butene-1 was absorbed in 
a solution of mercuric nitrate cata- 
lyzed by silver nitrate ‘?), the residual 
volume of n-butane was measured and 
the volume of butene-1 obtained by 
subtraction. 


Analysis for acetone 

Acetone reacts almost quantita- 
tively with hydroxylamine hydro- 
chloride in aqueous solution releasing 
hydrochloric acid, which may be 
titrated with standard sodium hydrox- 
ide solution. The size of the sample 
was chosen to give a reasonable volume 
of hydrocarbon gas (at least 30 ml.), 
so the amount of acetone present 
varied over a wide range. Normally a 
2% hydroxylamine hydrochloride 
solution was used. The volume of 
hydroxylamine hydrochloride solution 
in the bubbler varied from 250 ml. to 
1000 ml. depending on the quantity of 
acetone present. At least four times 
the theoretical amount of hydrochlor- 
ide required to react with the acetone 
was used in order to ensure that the 
reaction was as complete as possible. 

After all the sample had _ been 
passed through the solution and the 
hydrocarbon gas had been collected 
in the buret, the solution was titrated 
potentiometrically to pH 3.75 with 
standard sodium hydroxide solution. A 
similar blank solution was titrated to 
pH 3.75, and a correction was applied 
to the original titre. 

In calculating the weight of acetone 
present a factor was used to allow 
for the fact that the reaction was not 
quite complete. In a series of nine 
tests using a known weight of acetone 
made up in sample bombs with 
hydrocarbon present, the factor was 
determined as 1.0284 with a standard 
deviation of 0.0050. 


Analysis for water 

The sample was dissolved in pyri- 
dine, and transferred to a titration 
vessel] which was closed to exclude 
atmospheric moisture. The sample 
was titrated with Karl Fisher reagent 
using the “deap-stop” direct titration 
method “, In nearly anhydrous media, 
acetone reacts with methanol in the 
Karl Fischer reagent to form a ketal 
with the release of water. To inhibit 
this reaction, the methanol content of 
the reagent was reduced to a low 
value, viz. 50 ml]. methanol in a solu- 
tion of 900 ml. pyridine, 133 gm. 
iodine, and 70 ml. sulphur dioxide, 
and the titration was carried out 
rapidly in dilute solution in pyridine. 
Allowance was made in the result for 
the small amount of water in the 
pyridine used to dilute the sample. 
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Figure 5—-Activity coefficient values for isobutane in acetone-water solvent con- 

taining 11.0 Wt. % water (C, Hydrocarbon-Free) showing effect of varying 

temperature. Solid curves are predicted by binary Margules equation. Dotted 
curves are predicted by binary van Laar equation. 
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Figure 6—Activity coefficient values for cis-butene-2 in acetone-water solvent 

containing 11.0 Wt. % water (C, Hydrocarbon-Free) showing effect of varying 

temperature. Solid curves are predicted by binary Margules equation. Dotted 
curves are predicted by binary van Laar equation. 
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Figure 7—Activity coefficient values for n-butane at 150°F. In acetone-water 

solvent showing effect of varying solvent composition (C, Hydrocarbon-Free ). 

Solid curves are predicted by binary Margules equation. Dotted curves are 
predicted by binary van Laar equation. 


a 


+——+_+—______——_ 


89 






+—_—+—— 


| 
ool $4 + . _ 


16.8 WT. % WATER IN SOLVENT 


40 


ACTIVITY COEFFICIENT OF BUTENE-!, x 


2.0 


oO 0.05 0.10 0.15 020 0.25 0.30 


MOLE FRACTION BUTENE -1, x, 


Figure 8—Activity coefficient values for butene-] at 150°F. In acetone-water 

solvent showing effect of varying solvent composition (C, Hydrocarbon-Free). 

Solid curves are predicted by binary Margules equation. Dotted curves are 
predicted by binary van Laar equation. 


carbons in these tables were calculated 
by Equation 2: 


Activity coefficients of pseudo- 
binary mixtures 

Experimental results of pseudo- 
binary systems of hydrocarbon and 
solvent are given in Tables 2, 3, and 4. zPyi 
The activity coefficients of the hydro- 





where 


(Pi —P) (vi - B,) 
y = eecevee 3 
log 21 2.303 RT (3) 


The value z, is a _ correction 
factor which takes into account the 
non-ideality of the vapor and the 
effect of pressure on the liquid (®). 
Hydrocarbon vapor pressures were 
obtained from correlations by Gris- 
wold and Brooks ‘), Values of liquid 
hydrocarbon specific volumes and 
virial coefficients were interpolated 
from tables given by Mertes and 
Colburn “), and Welty, Gerster, and 
Colburn ‘), 


The experimental equilibrium data 
are shown in Figures 3 to 8, inclusive, 
as the logarithm of the hy drocarbon 
activity coefficient plotted against the 
mole fraction of hydrocarbon in the 
liquid phase. Figures 3, 4, 5, and 6 
show the effect of temperature on 
each hydrocarbon-solvent system, 
while Figures 7 and 8 show the effect 
of varying solvent composition at 
constant temperature for n-butane 
and butene-1, respectively. The solid 
curves through the experimental data 
indicated as circles and triangles were 
drawn using the binary Margules 
three-suffix equations (5), These may 
be written as: 


logy: = x2? [Ar + 2(Ae — Ai)xi]. .(4) 
limit log yi: asx1 ~O = Ai. 

1? [Az + 2(Ai — Agz)xe]. .(5) 

limit log y2 as x2 ~ O = 


log y2 = x 


The dotted curves through the ex- 
perimental data were drawn using the 
binary van Laar two-suffix equations 
(6) which may be written as: 


Ay * , 
w= = X2 
i." 


limit log y: as x1 ~ O = Ai. 


log yi = ae A 


logy2 = __ Aa xi? (7) 


Ao P 
E + . «| 


limit log y2 as x2 ~O = Ag 


In equations 4 to 7 inclusive, the 
subscript 1 refers to the hydrocarbon 
component, and the subscript 2 refers 
to the solvent mixture, acetone-water. 


In order to draw curves through 
the experimental points, values of the 
terminal activity coefficients, A, and 
A., must be known for the Margules 
and van Laar equations. The pro- 
cedure used to obtain values of A, 
and A, from the experimental activity 
coefficients was to apply the method 
of least squares to Equations 4 and 6. 
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TABLE 3 


PsEupo-Binary VApor-LiQuID EQUILIBRIUM DATA FOR n-BUTANE, BUTENE-1, ISOBUTANE, AND CIs-BUTENE-2 IN 
ACETONE-WATER SOLVENT CONTAINING 11.0 Wr. % WATER (Cy HyDROCARBON-FREE) 


Vapor Analysis (Mole %) Liquid Analysis (Mole %) Activity 
Total | ; a aes ___| Correction | Coefficient 


Run No. Pressure | Factor of 
Mm. Hg. | Acetone Water Hydrocarbon | Acetone Water Hydrocarbon z Hydrocarbon 
| 


en — Mian Ne ee 


n-Butane, Temperature 99.7°F. 


.083 10.54 
.083 an 
.083 .29 
.065 48 
.065 .90 
.064 soo 
.049 36 
.032 85 
.032 75 
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.022 .98 


mn 
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n-Butane, Temperature 149.7°F. 
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101 
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.037 
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.16 
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tm 
79. 
81. 
| 80. 


33 1474 |} 51. 
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| 
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Butene-1, Temperature 99.7°F. 
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Butene-1, Temperature 149.7°F, 
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Butene-1, Temperature 199.0°F. 
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n-Butane, Temperature 199.0°F. 
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TABLE 3 (Continued) 


PsEubDO-BINARY VAPOR-LIQUID EQUILIBRIUM DATA FOR n-BUTANE, BUTENE-1, ISOBUTANE, AND CIs-BUTENE-2 IN 
ACETONE-WATER SOLVENT CONTAINING 11.0 Wr. % WATER (Cy HYDROCARBON-FREE) 
Vapor Analysis (Mole %) Liquid Analysis (Mole %) Activity 
Total _ ______} Correction Coefficient 
Run No. | Pressure Factor of 
Mm. Hg Acetone Water Hydrocarbon | Acetone Water Hydrocarbon Zz Hydrocarbon 
Isobutane, Temperature 99.7°F. 
62 1257 24.1 2.9 73.0 ooo «| 2h2 2.83 1101 | 9.64 
1259 22.1 2.9 75.0 69.5 | 27.6 2.86 1.101 | 9.81 
61 2031 15.4 ey 82.9 67.0 26.8 6.22 1.069 7.82 
2028 14.4 Lo 83.7 66.8 26.9 6.27 1.069 7.81 
60 2527 11.6 1.4 87.0 64.7 29.9 9.83 1.048 6.33 
2523 11.0 1.8 87.2 64.2 23.4 10.06 1.049 6.19 
59 2797 9.0 1.8 89.3 61.7 25.6 12.69 1.038 5.51 
2780 10.2 2.0 87.8 62.2 24.9 12.94 1.038 5.28 
Isobutane, Temperature 149.7°F. 
87 1793 | 44.7 | a4 47.6 70.0 28.4 1.59 1.184 8.62 
1796 42.9 7.9 49.2 69.7 28.7 1.62 1.184 8.76 
86 2767 23.4 5.6 71.0 68.3 22 3.96 1.150 7.74 
2764 23.8 5.6 70.6 68.4 27.7 3.90 1.150 7.81 
88 3481 20.0 4.2 75.8 66.8 26.9 6.31 1.126 6.39 
85 4288 14.2 oz 82.6 | 64.9 25.9 9.18 1.098 o.75 
cis-Butene-2, Temperature 99.7°F. 
95 744 39.6 5.1 55.9 67.8 28.7 o.04 1.066 5.33 
744 39.4 4.9 55.7 67.8 28.7 3.48 1.066 5.42 
94 1155 23.9 5.2 72.6 64.9 26.6 8.56 1.049 4.39 
1152 24.. 3.9 4258 65.1 26.2 8.66 1.049 4.30 
96 1471 iv 2.8 79.4 61.3 24.2 14.52 1.036 3.90 
1476 18.0 3.0 79.0 61.0 24.7 14.34 1.036 Sian 
118 1611 17.5 $3.2 79.3 58.4 23.9 17.68 1.030 317 
1611 17.0 2.9 80.1 58.0 24.1 17.93 1.030 3.16 
cis-Butene-2, Temperature 150.0°F. 

117* 1506 50.9 9.3 39.8 | 70.1 nt so 2.57 | <P eae of a. 87 
1506 51.6 8.9 39.5 | 68.9 28.5 2.60 1.116 | 5.06 
91 2273 29.7 6.0 64.3 | 65.5 21.3 7.40 1.089 4.29 
| 2265 27.8 6.2 66.0 66.2 26.3 7.47 1.089 4.35 
92 | 3004 20.4 4.5 75.1 61.5 24.6 13.92 1.064 3.44 
| 2986 21.3 5.0 73.7 61.6 24.2 14.16 1.065 3.30 
98 3285 16.8 4.5 78.7 59.7 22.1 18.25 1.055 2.99 
3285 7.5 oo 79.0 59.8 ene 18.08 1.055 3.03 


*The temperature for Run 117 was 150.4°F. 





The terminal activity coefficients are 
tabulated in Table 5. 

A pseudo-component — treatment 
may be used when the ratio x,, 
Xwater (K) remains constant. The 
analytical results in Tables 2 to 4 
show some variation in the value of 
K. The resulting variation in hydro- 
carbon activity coefficients is small, 
and the values of the terminal activity 
coefficients are essentially those for a 
solvent with the average value of K. 
An alternative treatment (4), in which 
hydrocarbon concentrations were 
estimated for an exact value of K, 
gave values up to 3% different for A 
and 6% different for A, in the worst 
case. The consistency of results, as 
shown in Figures 3 to 8 inclusive, 
indicates that the C, hydrocarbon- 
acetone-water pseudo- -binary systems 
can be reliably interpreted by the 
Margules and van Laar equations. 
Further evidence will show that treat- 


Pn 
etone/ 


26 


ment of acetone-water as one com- 
ponent in a pseudo-ternary system 
with two hydrocarbons can also be 
interpreted by Margules and van Laar 
equations. 


Calculation of enthalpy values 

For design and optimum operation 
of extractive distillation towers know]- 
edge is required of heats of solution 
for calculation of enthalpies and fluid 
rates within the tower. Since the ex- 
perimental data on activity coeffi- 
cients appear consistent, reliable values 
of heats of solution may be derived 
using Equation 8. 
Li 


= .(8) 
2.303R 


d log y; 
d(1/T) 


The function on the left side of 
Equation 8 is the slope of a curve 
formed by plotting the logarithm of 
the activity coefficient for a given 
hy drocarbon of composition x, against 


the reciprocal of the absolute tem- 
perature. When x, = 0, this involves 
plotting the terminal activity coefh- 
cient A, against 1/T. The A, values 
of Table 5 appear to be linear func- 
tions of 1/T, from which it follows 
that the heat of solution for the four 
hydrocarbons in the solvent is essen- 
tially independent of temperature in 
the range of 100° to 200°F. 

While A, can be calculated accur- 
ately by applying experimental results 
to Equations 4 and 6, the determina- 
tion of A, is much less accurate, 
particularly at high temperatures 
where the range of hydrocarbon con- 
centrations was limited. - alues of A, 
were plotted against 1/T. For n-bu- 
tane, isobutane and aa straight 
lines were drawn through the points 
at 99.7° and 149.7°F., while for cis- 
butene-2 a line parallel to that for 
butene-1 was drawn through the 
point at 99.7°F. Smoothed Margules 
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TABLE 4 


PsEupo-BINARY VAPOR-LIQUID EQUILIBRIUM DATA FOR n-BUTANE AND BuTENE-1 IN ACETONE-WATER SOLVENT 
CONTAINING 16.8 Wr. % WATER (Cy HyYDROCARBON-FREE) 





Vapor Analysis (Mole %) Liquid Analysis (Mole %) Correction Activity 
Total | i —_ . - : Factor Coefficient 
tun No. Pressure | =————<CStsi<CSsté‘ Zz of 
— a. the | Acetone Water Hydrocarbon Acetone Water Hydrocarbon Hydrocarbon 
See silat a eel 
n-Butane, Temperature 150.3°F. 

x tt7z 0 CO | 483 8.4 49.3 59.7 38.8 1.43 1.128 12.58 
” 1772 | 41.3 8.5 50.2 59.7 38.9 1.41 1.128 1299 
106 2535 26.7 5.9 67.4 59.2 37.7 3.08 1.101 11.14 

2535 26.9 5.9 | 67.2 59.4 37.5 3.06 1.101 11.19 
104 3037 19.3 5.0 75.7 58.3 37.3 4.37 1.083 10.40 
3036 19.2 5.0 75.8 57.4 38.0 4.53 1.083 10.04 
101 3466 17.7 3.4 78.9 56.4 37.4 6.21 1.068 8.59 
3466 | 17.4 wl 79.5 56.4 37.4 6.26 1.068 8.59 
Butene-1, Temperature 150.4°F. 
1507 | 41.5 9.5 48.9 59.7 38.8 1.51 1.160 8.70 
63 1507 | 43.3 8.8 47.9 59.6 38.8 1.63 1.160 7.88 
115 2021 | 34.7 73 58.1 58.3 38.9 2.76 1.142 746 

2021 33.6 7.4 59.0 58.8 38.4 2.79 1.142 7.50 

112 2519 21 =| «6.9 68.0 57.9 37.8 4.34 $,925 6.83 
2519 26.1 | 5.7 68.1 58.1 37.5 4.41 1.125 6.73 
109 3028 19.6 5.6 74.8 56.9 36.5 6.59 1.109 5.86 


values of A, and A, are presented in 
——T Table 6. The linear relationships of 
A, and A, with 1/T was chosen to 
simplify interpolation of values at 
intermediate temperatures. Such a 


TERMINAL ACTIVITY COEFFICIENTS OF PSEUDO-BINARY SYSTEMS 











Terminal Activity Coefficients relationship is consistent with data on 
Average aaa — — C, hydrocarbon-furfural systems (1.4), 
oon Temp. | Margules van Laar and there is no reason to believe a 
Dinary Synem Seiveat, oF ee Ye more complicated relationship exists 
: wt. % in the present system. Nevertheless, 
. Ay ; As Ay : fs = caution should be applied when using 
n-Butane-Solvent 10.8 99.7 1.038 0.677 1.039 0.755 values of the constants at 200°F. 
5.4 150.3 0. 775 0.709 0. 773 | 0.718 The binarv Margules equation, 
11.0 149.7 | 0.978 | 0.603 | 0.978 | 0.692 Ecuati Ape. 7 
16.8 150.3 1.163 | 0.400 | 1.166 | 0.676 uation 4, with smoothed values of 
10.9 199.0 0.890 | 0.492 | 0.891 | 0.603 A, and A, from Table 6 was used to 
| aoe se r compute activity coefficients at var- 
Isobutane-Solvent 0 99.7 | 1.070 Hee aaa on ious hydrocarbon concentrations. 
, or oa | These values plotted against 1/T 
Butene-1-Solvent | 10.9 | 99.7 | 0.796 | 0.695 0.793 | 0.715 gave straight lines for each hydro- 
| 5.4 | 150.4 | 0.585 0.623 0.585 | pape carbon concentration. The slopes of 
| 10.9 | as ae ge car ec the lines were determined, and by 
a | 199.0 0.725 0.398 0.721 | 0.500 means of Equation 8 the differential 
. heats of solution of the liquid hydro- 
cis-Butene-2-Solvent 11.2 99.7 0.791 0.650 0.791 0.667 carbons in the acetone-water solvent 
10.9 150.0 0.741 eal oa One were calculated. These are presented 
in Figure 9. The curves in Figure 9 
were integrated graphically to obtain 
integral heats of solution, which are 
ai presented in Table 7. The sign is 
TABLE 6 positive since these heats of solution 
SMOOTHED MARGULES PsSEUDO-BINARY TERMINAL ACTIVITY COEFFICIENTS IN represent the heat absorbed during 
SOLVENT CONTAINING 11.0 Wr. % WATER (Cy HyDROCARBON-FREE) solution of liquid hydrocarbon in the 
solvent. 
n-Butane Isobutane Butene-1 cis-Butene-2 Similar results for heats of solution 
ON eset ns wis wo) pi ete would be obtained using Equation 6 
- | _ - - \ \ fe ie he and smoothed values of the van Laar 
sccancistour =e eres icen a : |" constants, A, and A,. As seen in 
99.7 | 1.047 | 0.677 | 1.070 | 0.642 | 0.796 | 0.695 | 0.791 | 0.650 Table 5 the A, values are a 
‘ i ic | > Margules values. 
149.7 | 0.964 | 0.603 | 0.986 | 0.598 | 0.756 | 0.610 | 0.747 | 0.565 — hfe = foe highen 
199.0 | 0.895 | 0.542 | 0.916 | 0.561 | 0.722 | 0.538 | 0.710 | 0.493 but essentially parallel to the Margules 


values when plotted against 1/T. 
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Figure 9—Differential heats of solution of liquid C, hydrocarbons in acetone- 
water solvent containing 11.0 Wt. % water (C, Hydrocarbon-Free) in temperature 


range of 100° to 200°F. 


TABLE 7 


INTEGRAL HEATS OF SOLUTION OF LiguID Cy HYDROCARBONS IN ACETONE-WATER 
SOLVENT CONTAINING 11.0 Wr. % WATER (Cy, HYDROCARBON-FREE) IN THE 
TEMPERATURE RANGE OF 100° To 200°F. 


Composition 
Mole % C, 
Hydrocarbon 


Integral Heats of Solution, Btu./Lb. Mole of Hydrocarbon Dissolved 





n-Butane Isobutane Butene-1 cis-Butene-2 
0.  2s70—Cté«‘(L;Ssé‘éilkwd‘SSC*dTSCt(‘«éaOSC*«*dSC<«*‘i«éi TIN 
5 2440 2430 1260 | 1360 
10 2300 2250 1260 | 1350 
15 2170 2090 | 1250 1330 
20 2050 1940 | 1240 | 1310 
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Since the values of the differential 
heats of solution are obtained from 
slopes of log y, versus 1/T, the results 
will be the same using both Margules 
and van Laar constants. Easier manip- 
ulations with the Margules equation 
was the reason for its use in these 
calculations. 


The saturated liquid enthalpy of a 
mixture is the sum of the enthalpies 
of the pure liquid components and the 
integral heat of solution of the liquid 
components when mixed together, 
and will be expressed at Bru. per 
pound of solution. A reference state 
was arbitrarily chosen as the saturated 
liquid at 0°F. Enthalpies of the pure 
liquid hydrocarbons were obtained 
from data published by Maxwell (19), 
while enthalpies of the pure liquid 
mixture of acetone and water were 
obtained from Shell Development 
Company “!!), These enthalpies were 
combined with the integral heats of 
solution to give saturated liquid 
enthalpy-concentration curves as 
shown in Figure 10 for the n-butane- 
acetone-water system. 


For multicomponent liquid or va- 
por mixtures, such as several C, 
hydrocarbons dissolved in acetone- 
water, enthalpies are conveniently ob- 
tained from the partial enthalpies of 
the components. The enthalpy is 
obtained by adding the products of 
weight fraction and partial enthalpy 
of all the components as indicated by 
Equation 9. 


h = Dwihi.............(9) 


Partial liquid enthalpies are deter- 
mined from saturated liquid enthalpy- 
concentration curves such as shown 
in Figure 10 by the method of inter- 
cepts (2), A tangent, for example 
A BC, is drawn to a curve at a given 
n-butane concentration, say 20 wt. %. 
Values of partial enthalpies are read 
at both extremities of the diagram. At 
C, the value is the partial enthalpy of 
n-butane at the given concentration 
of 20 wt. % and a temperature of 
150°F. At A, the value is the partial 
enthalpy of the solvent at 150°F. con- 
taining 20 wt. % n-butane. Values 
such as C yield partial liquid enthalpy- 
concentration curves for the four 
hydrocarbons studied. These curves 
are presented in Figure 11. Values 
such as A are found to be almost 
identical for all four hydrocarbon- 
solvent systems at given hydrocarbon 
concentrations and temperatures. 
These values are averaged and pre- 
sented in Figure 12 as partial liquid 
enthalpy-concentration curves for 
acetone-water solvent containing any 
C, hydrocarbon. 


The partial enthalpies of vapor 
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Figure 10—Saturated liquid enthalpy- 

concentration curves for the n-butane- 

acetone-water system. Solvent contains 

11.0 Wt. % water (C, Hydrocarbon- 
Free). 


phase mixtures of C, hydrocarbon- 
acetone-water are the same as the 
enthalpies of the pure components at 
any given temperature and solvent 
composition. This result follows from 
the assumption that the heat of mixing 
of a C, hydrocarbon and solvent in 
the vapor phase is negligible. 


Activity coefficients of pseudo- 
ternary systems 

The pseudo-ternary system of n- 
butane-butene-1-acetone-water was 
investigated at 100° and 150°F., at a 
water concentration of 11.0 wt. % 
(C, hydrocarbon-free). Values of 
activity coefficients for n-butane and 
butene-1 were computed by Equa- 
tions 2 and 3. Results are given in 
Table 8. 

An analysis of theoretical equations 
which describe ternary vapor-liquid 
equilibrium data has been made by 
Wohl (13), The ternary Margules 
three-suffix equations, for example 
Equation 10 for component 1, and 
the ternary van Laar two-suffix equa- 
tion, for example Equation 11 for 
component 1, were used in the present 
investigation. Similar equations may 
be written for the other components 
by the cyclic rearrangement of sub- 
scripts, that is, substituting 2 for 1, 
3 for 2, and 1 for 3. 


log vi = X2? [Ais + 2x; (Ao-1 — Aj-2)] + 


ms" [Ais + 2x1 (A3-1 * A1-3)] + 
X2X3 [Ae-1 + Ai-3 — A3-2 + 
2x1 (A3-1 = A1-3) + 2x3 (A3-2 = Ao-3) = 


Cyo3 (1 — 2x1)]. .(10) 
(see Equation 11 on page 30) 


For the pseudo-ternary system 
studied, the subscripts 1, 2, and 3 refer 
to n-butane, butene-1, and acetone- 
water solvent, respectively. The ter- 
minal activity coefficients, A,.. and 
A,.,» for the n-butane-butene-1 sys- 
tem were assumed to be zero, assum- 
ing that a binary system of two 
hydrocarbons is ideal ‘?). In the van 
Laar equation, where these constants 
appear as a quotient, the value may 
be computed in terms of the other 
four binary constants of the ternary 
system (6), Equa-:ion 12 expresses the 
relationship. 





Figure 11—Partial liquid enthalpy of 

C, hydrocarbons in acetone-water sol- 

vent containing 11.0 Wt. % (C: Hydro- 
carbon-Free ). 


BTU/LB. SOLUTION. 


TOTAL ENTHALPY, 


PARTIAL ENTHALPY OF HYDROCARBON, BTU/LB. OF HYDROCARBON IN SOLUTION. 





200 °F 


1SO°F A 


; 





1\00°F 





90 100 
AMOUNT N-BUTANE DISSOLVED IN SOLVENT, WT.%. 


180 





160 

200°F 
140 
120 

ISO°F 
100 
™ 10°F 


—4— - N-BUTANE 

60 -O— - ISOBUTANE 
-~O—-- BUTENE-! 
~O—-- CIS-BUTENE-2 





40 at He 
0 5 10 15 20 25 30 
AMOUNT G4 HYDROCARBON DISSOLVED IN SOLVENT, 
WT.%. 


The Canadian Journal of Chemical Engineering, February, 1958 29 








TABLE 8 


PsEuDO-TERNARY VAPOR-LIQUID EQUILIBRIUM DATA FOR THE SYSTEM n-BUTANE-BUTENE-1-ACETONE-WATER 
SOLVENT CONTAINING 11.0 Wr. % WaTER (Cy HYDROCARBON-FREE) 


Potal 


Pressure 
Mm. Hg. 


2007 
2001 
1999 
2009 
2012 
2755 
2762 
2765 
2782 
3515 
3525 
3501 
3547 
4041 
4041 
4068 
4073 
4078 


*The temperature for 


Component 


cis-Butene-2 


Isobutane 


Equation 11 


Reo es w 


So) 


Sonam os 


= 


Vapor Analysis (Mole %) 


Acetone Water n-Butane 


.0 Fe 19.6 
a 5.4 16.5 
S 2.4 18.7 
8 4.1 29.9 
3 3.7 29.5 
8 ace 57.6 
Zi ont 
2.5 ze.¥ 
es 62.0 
0 2.3 24.3 
Za 24.3 
eis 48.5 
mae 63.4 
7.2 17.4 
6.8 16.9 
6.8 16.5 
7.4 38.8 
can 38.5 
5.9 a2 
1 ee 
5.3 22.3 
5.1 cs.7 
8 3.9 23.4 
2 3 23.3 
8 a7 22.4 
7 wee 52.8 
3.1 $4.3 
2.4 54.1 
3.0 20.4 
3.0 27.3 
3.0 26.3 


Runs 119 and 124 was 15 


SAMPLE CALCULATION O1 


Liquid 
Phase 
Composition 
xX, Mole 
Fraction 
0.0195 
0.0331 
0.0174 
0.0258 


0.9042 





Te Ww 


mUuUmumnNnNNuUuN 
CO = 


mmw 


Acetone 


70. 
68... 
66. 
65. 
66. 
06. 
63. 
63. 


63 


60. 
59. 
59. 
59. 


70. 
67. 
67. 
67. 
66. 
62. 


\mount of 
Component 
Relative to 
Butene-1, ©; 


100.0 


a 


NON EA Ws 


CoOnT Iw 


No DW Wwe ee 


=~ 


Liquid Analysis (Mole %) 


Temperature = 99.7°F. 


0. 


‘emperature = 149.7°F. 


Water 

anak 

27.4 0 

25.6 

26.8 

26.5 

26.6 

a3.7 

23.4 

25.2 

i 

22.0 

eect 9. 

Zi.S 1 

28.1 0. 

27.8 

23.4 

26.9 

26.7 

ee | 

25.6 

25.3 

ete 

26.3 

26.0 

26.4 

26.0 

24.3 

24.3 ¢ 

24.1 4. 

23.9 4. 

23.9 4, 
TABLE 9 


VAPOR PHASE COMPOSITION BY THE I 


Volatility 
of Component 
Relative to 
Butene-1 at 
150.4°F. 


a 
1.000 
0.778 
1.275 


1.695 


0.0353 


n-Butane 


ees we er 


Nm 


0. 
0. 


Om WW We ee N 


oo 


Activity Activity 
Correction | Coefficient | Correction | Coefficient 
Factor of Factor of 
Butene-1 Zi n-Butane Z2 Butene-1 
1.61 1.083 9.85 1.097 O.ee 
3.55 1.070 8.44 1.085 S43 
6.36 1.052 1.89 1.069 4.56 
4.97 1.052 7.50 1.069 4.57 
2.22 1.051 7a 1.068 4.51 
2.25 1.051 7.52 1.068 4.65 
9.95 1.034 5.42 1.052 3.69 
9.95 1.034 5.44 1.052 3.70 
3.48 1.030 5.91 1.049 3.95 
14.39 1.021 4.72 1.040 2.99 
14.12 1.022 4.76 1.041 3.01 
8.60 1.019 4.39 1.039 ouae 
5.89 1.019 4.15 1.038 2.87 
2.67 1.119 7.64 1.141 5.20 
2.61 1.119 8.10 1.141 4.99 
2.66 1.119 7.61 1.142 4.95 
1.19 1.119 cae 1.141 5.10 
1.21 1.119 1.4% 1.141 5.02 
5.11 1.093 6.70 1.118 4.31 
5.13 1.093 6.83 1.118 4.43 
a. 1.093 6.51 1.118 4.37 
2.81 1.092 6.96 1.117 4.47 
8.51 1.065 5.28 1.091 3.83 
8.59 1.064 5.26 1.091 3.68 
8.58 1.065 5.38 1.092 3.71 
3.87 1.064 5.38 1.090 3.76 
5.96 1.047 4.51 1.074 3.12 
5.85 1.047 4.37 1.074 3.07 
1.68 1.046 4.35 1.073 3.09 
1.63 1.046 4.57 1.073 3.20 
1.78 1.046 4.25 1.073 $.2i 
RELATIVE VOLATILITY METHOD 
Calculated Experimental 
Equilibrium Equilibrium 
Vapor Phase Vapor Phase 
ax Composition, Composition, 
y, Mole vy, Mole 
Fraction Fraction 
0.0195 0.136 0.136 
| | 
0.0258 | 0.180 | 0.185 
0.0222 | 0.155 | 0.151 
0.0437 | 0.305 | 0.305 
| | 
0.0319 0.223 | 0.223 
| 
0.1431 | | 


At a first approximation the ternary 
constant, C,,,, in the Margules equa- 
tion was assumed to be zero. As will 
be seen later, this assumption was 
justified. 

In Figures 13, 14, 15 and 16 the 
experimental values of hydrocarbon 


(88) activity coefficients, indicated by 


circles, are shown as functions of 
composition. The abscissa is mole 
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ACTIVITY COEFFICIENT OF N-BUTANE , 





fraction of hydrocarbon component, 
100 x; on 100 xs 


X2+X3 Xi +Xs 
represents the relative —— of the 
other two components. Values of 
activity coefficients predicted by 
Equations 10 and 11 are shown as 
solid and dotted curves, respectively. 
Smoothed values of ‘the terminal 
activity coefficients were used. Good 
agreement was obtained between ex- 
perimental activity coefficients and 
those predicted by the ternary equa- 
tions, as shown in Figures 13 to 16, 
inclusive. Hence, a value of C,,, =O 
appears to be justified. 

It is not surprising to see both the 


van Laar and Margules ternary equa- 


while the parameter 

















tions describe equally well the present 
system at 100 and 150 F. In systems 
like the present or the C, yy drocarbon- 





furfural system (1,3,4) which are 80 ices Se 
slightly dissimilar as indicaved by the 
10.0 | ltOoo°F 
98.4 | o 
9.0 a 60F + - T = Pate T 
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6.0 S Figure 12—Partial liquid enthalpy of acetone-water solvent containing 11.0 Wt. 
\ 96.2 % water (C, Hydrocarbon-Free) in solution with any C, hydrocarbon. 
oN 

\ 897 
oor Nae | | ; 


ratio of A, to A, being near 1 both 
equations give similar results (*), Thus, 
calculations involving the ternary 
equations justified using the Margules 
equations because they were easier to 
work with. Since the ternary equations 
were not checked experimentally at 
200°F. caution should be exercised 
when using the terminal activity co- 
efficients at this temperature. 


N. 











Relative volatility values 


In distillation applic ations relative 
volatility is of prime importance. The 
dependence of relative volatility on 
other variables is often expressed ir 
terms of activity coefficients as pes 
fined by Equation 13. 


ACTIVITY COEFFICIENT OF N- BUTANE 


P\Z2 
ai-2 = yu =) «st RS) 
¥2 Pozi 


Figure 13—Activity coefficients of n- 
butane in mixtures with butene-l and 
acetone-water solvent at 100°F. Sol- 
vent contains 11.0 Wt. % water (Cy 
Hvydrocarbon-Free). Solid curves are 
predicted by ternary Margules equa- 
tion. Dotted curves are predicted by 
° 0.05 0.10 0.15 0.20 0.25 0.30 ternary van Laar equation. Values of 
MOLE FRACTION N- BUTANE, X,. parameter are 100 x;/ (x2 + xs). 
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The ratio y,;/y2: may be calculated 
fron) the ternary Margules or van 
Laar equations. The former equations 
yield the simpler expression, given by 
Equation 14, assuming that A,_,, 


A,.,, and C,,, are each zero. 


log = = x3 [(1 — 2x1) Ai-3 — X3Ae-2 + 


2x1A3-1 — (1 ~ Xs) A3-2]. (14) 


Values of vapor pressures were ob- 
tained from Griswold and Brooks “), 
while the correction factors z,, Z., 
were calculated by Equation 3. 
Smoothed values of the terminal 
activity coefficients were used in 
Equation 14. 

Volatilities of n-butane isobutane, 
and cis-butene-2 relative to butene-1 
present in solvent containing 11.0 wt. 
% water (C, hydrocarbon-free), cal- 
culated by Equations 13 and 14, are 
shown in Figure 17, 18 and 19 as 
functions of temperature and compo- 
sition. Figure 20 shows the volatility 
of the solvent mixture containing 11.0 
wt. % water (C, hydrocarbon-free) 
relative to butene-1 as a function of 
temperature and composition. In 
Figures 17, 18 and 19 the relative 
volatility is seen to decrease with 
increasing temperature and with in- 
creasing total hydrocarbon concentra- 
tion in the liquid phase. 

In a multicomponent extractive 
distillation calculation the composi- 
tion of vapor phase in equilibrium 
with a liquid phase is conveniently 
determined by the relative volatility 
method “'#), Equation 15 is used in 
the method. 


Qn Xp 
Wig Se ena weno (15) 


Lax 


It is desirable to have the volatilities 
relative to one of the components; 
butene-1 was chosen in this investiga- 
tion. 

Table 9 provides a sample calcula- 
tion of the vapor-phase composition 
by the above method when the liquid 
phase composition is known. The 
liquid phase composition was ob- 
tained from a run at 150.4°F. with a 
multicomponent C, hydrocarbon and 
solvent mixture containing 11.0 wt. % 
‘vater (C, hydrocarbon-free). The 
experimental vapor phase composition 
is also shown to compare with the 
calculated composition. The agree- 
ment is satisfactory and indicates that 
the relative volatility data of Figures 
17 to 20 inclusive, may be used in 
calculations provided some caution is 
exercised in the 200°F. range where 
the bulk of the smoothing has been 
done . 
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Figure 16—Activity coefficients of butene-l in mixtures with n-butane and 
acetone-water solvent at 150°F. Solvent contains 11.0 Wt. % water (C, Hydro- 
carbon-Free). Solid curves are predicted by ternary Margules equation. Dotted 
curves are predicted by ternary van Laar equation. Values of parameter are 


100 x;/(x: 


Conclusions 

Vapor-liquid equilibrium data for 
C, hydrocarbons-acetone-water sys- 
tems have been measured. The con- 
sistency of the experimental data 
interpreted as activity coefficients 
confirm the reliability of the flow- 
type apparatus and the analytical 
methods. Correlation of the experi- 
mental pseudo-binary and _ pseudo- 
ternary data by use of binary and 
ternary Margules three-suffix and van 
Laar two-suffix equations was satis- 
factory in the range of composition 
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T Xs). 


and temperature studied. The solvent, 
acetone-water of constant composi- 
tion, can be treated as one component. 
Treatment of the vapor-liquid equi- 
librium data with Margules three- 
suffix equations to yield enthalpy and 
relative volatility data was found to 
be adequate. The enthalpy and rela- 
tive volatility data are expressed as 
functions of temperature and com- 
position. Such data can be used for 
calculation of rates of flow and com- 
positions of vapor and liquid in 
extractive distillation towers. 
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Nomenclature P 
Ay = terminal activity coefficient of Pp 
hydrocarbon in binary equa- ; 
tions, limit log y:, as x; — O. 
. Og aig ) 
Ao = terminal activity coefficient of I 
solvent in binary equations 
limit log y2 as x2 > O. R 
\i-2 etc. = terminal activity coefficient in 
ternary equations, limit log 
Yiasx; > Oand x: > 1. , 
B = second virial coefficient in the 
equation of state, ml./gm.mole. V 
i. = compressibility factor. 
Cyo3 = ternary constant in Margules Vv 
three-suffix equation for ter- 
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nary system of 
lL, 2, ana Ss. 


components 


= differential operator. 


total enthalpy of liquid, Btu./ 
lb. of solution. 


partial enthalpy of compon- 
ent in solution, Btu./lb. of 
component. 


ratio of mole fraction acetone/ 
water in liquid phase. 

differential molal heat of solu- 
tion of liquid Cs hydrocarbon 
in acetone-water solvent, 
Btu./lb. mole of hydrocarbon. 


= molecular weight. 
= weight of hydrocarbon, gm. 
= total pressure of system, mm. 


Hg. 

vapor pressure of pure com- 
ponent, mm. Hg. 

partial pressure of compon- 
ent, mm. Hg. 

gas law constant, 62,365 
(ml.) (mm. Hg.) (°K.) (gm. 
mole) or 1.987 Btu./(°R.) (Ib. 
mole). 

absolute temperature, °R. or 
"i. 

volume of hydrocarbon gas, 
ml. 

liquid specific volume, ml./ 
gm. mole. 


—$$____}_— 


140 
TEMPERATURE, °F 





w = weight fraction of component, 

x = mole fraction of component 
in liquid phase. 

Vv = mole fraction of component in 
vapor phase. 

z = gas law deviation and liquid 
pressure factor of component. 

ai-2 = relative volatility of compon- 
ent 1 to component 2. 

Y = activity coefficient in liquid 
phase. 

> Ronee 

2 = summation sign. 

Subscripts 

1 = component 1, any Cy, hydro- 
carbon in binary equations or 
n-butane in ternary equations. 

2 = component 2, solvent in bin- 
ary equations or butene-1 in 
ternary equations. 

3 = component 3, solvent in ter- 
nary equations, 

i = components 1 to i. 

n = any component. 
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Figure 17—Relative volatility of n-butane to butene-1 in acetone-water solvent as a function of temperature and composi- 


tion. Solvent contains 11.0 Wt. % 


water (C, Hydrocarbon-Free). Parameter is total hydrocarbon dissolved, Mole %, in the 


liquid phase. Dashed, solid, and dotted curves are for 0, 50, and 100 Mole % n-butane, respectively, on a solvent-free basis. 
Curves are predicted by Margules ternary equations. 
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Figure 19—Relative volatility of cis-butene-2 to butene-1, in acetone-water solvent as a function of temperature and 

composition. Solvent contains 11.0 Wt. % water (C, Hydro-carbon-Free). Parameter is total hydrocarbon dissolved, mole 

%, in the liquid phase, dashed, solid, and dotted curves are for 0, 50, and 100 mole % cis-butene-2, respectively, on a 
solvent-free basis. Curves are predicted by Margules ternary equations. 
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Figure 20—Relative volatility of acetone-water solvent to butene-1 as a function of temperature and composition. Solvent 


contains 11.0 Wt. % 
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The Leaching of Manganese 


from Pyrolusite Ore by Pyrite’ 


A new method of leaching low-grade manganese ore (< 5% 


G. THOMAS, 


and B. J. P. 


Mn) 


utilizing pyrite as a reagent was developed. By heating an aerated slurry 
of ore and pyrite and by recycling the leach liquor at controlled acidity 
to attack fresh ore, a manganous sulphate solution free of iron was 
obtained which when heated to about 175°C., yielded manganous sulphate 
monohydrate. During the leaching process, iron in the liquor was alter- 
nately oxidized by the pyrolusite and reduced by the pyrite. In an 
alternative leaching procedure, the heating of an aqueous slurry of 
pyrite in an autoclave first under oxygen pressure and then in the 
absence of oxygen produced an acidic ferrous sulphate solution suitable 
for the subsequent leaching of pyrolusite ore at ambient temperatures. 


I conjunction with a major iron 
mining operation in Ontario, a 
large annual tonnage of manganiferous 
low-grade iron ore is brought to the 
surface as waste material. The average 
composition of this material is 50% 
silica, 40% iron oxides, 6% clay and 
accessory minerals and 4% manganese 
oxides, principally pyrolusite, MnO, 
If the manganese could be recov ered 
from this waste material, the amount 
available would almost satisfy Canada’s 
annual manganese requirement of 
30,000 tons. During a Mines Branch 
investigation of this material, a new 
process for leaching manganese was 
developed. This paper describes the 
development and chemistry of the 
leaching process. 

There is an abundant literature on 
processes for the recovery of manga- 
nese from relatively high- -grade pyro- 
lusite ore. Of the many processes, the 
most interesting approach for the re- 
covery of manganese from this low- 
grade material a»peared to be that 
based on the use of pickle liquor from 
the steel industry. Hoak and Coull 


have shown that pickle liquor, which 
contains free sulphuric acid and fer- 
rous sulphate, would give a 98% 
extraction from ores containing 15 to 
25% manganese. There are many 
literature references (7-5) stating that 
sulphuric acid and ferrous sulphate 
are formed, in the presence of water, 

by the air oxidation of the relatively 
cheap and abundant mineral pyrite, 
FeS,. It seemed reasonable therefore 
to determine if manganese could be 
leached from an aerated slurry of 
pyrolusite ore using pyrite as the 
leaching agent. 


Experimental leaching of 
manganese dioxide 


The initial experimental work to 
test the feasibility of leaching pyro- 
lusite ore by pyrite was done using 
pure manganese dioxide. When a 
slurry of pyrite (80% minus 200 
mesh) and manganese dioxide in 
water was treated at 175°C. and under 
100 psig oxygen pressure for two 
hours in an autoclave, very little 
manganese was dissolved. Although 


1Manuscript received December 9, 1957. 
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manganese dioxide is insoluble in 
dilute sulphuric acid, the addition of 
this acid was found to accelerate the 
dissolution of manganese dioxide by 
pyrite. For example, when 1% v/v 
H.SO, was substituted for water in 
the above experiment, the manganese 
was dissolved completely in two 
hours. 


Additional experiments were done 
to determine the effect of lower tem- 
peratures on the dissolution of man- 
ganese dioxide. The experimental 
conditions and results for these ex- 
periments, given in Figure 1, curve 1, 
show that the dissolution of manga- 
nese by pyrite in 1% v/v H SO, was 
complete after two hours at 110° eS 
Analysis of the leach liquors from 
these experiments showed that high 
manganese:iron ratios of from 8:1 to 
42:1 were obtained. 

A pachuca was used for many of 
the subsequent experiments. It con- 
sisted of a glass cylinder fitted with a 
sintered glass disc in its conical base 
for air dispersion. In an experiment 
using the pachuca, air at atmospheric 
pressure was bubbled through a 
slurry of 2 gm manganese dioxide 
and 2 gm pyrite in 500 ml 2% v/v 
H,SO, at 90°C. Under these leaching 
conditions with the open vessel, the 
MnO, was completely dissolved in 
six hours. 


Experimental leaching of 
pyrolusite ore 


To determine if pyrite could be 
used to leach manganese from the 
low-grade pyrolusite ore, samples of 
ore containing 3.7% Mn were treated 
at various temperatures in an auto- 
clave for two hours with 1% v/v 
H,SO, and 100 psig oxygen pressure. 
The resales of these experiments are 
shown in Figure 1, curve 2. A com- 
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Figure 2—Rate of extraction of manganese from ore in the presence of pyrite. 


parison of these results with those for 
the leaching of manganese from pure 
MnO., curve 1, shows that the leach- 
ing of manganese from the ore was 
somewhat slower but nevertheless 
promising. 
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Additional experiments on_ the 
leaching of pyrolusite ore were then 
done in an open vessel. The leaching 
experiments were done with pyrite 
and 2% v/v H,SO, in a pachuca at 
temperatures of 80 C and 90°C. and 


with acid to ore ratios of 3:1 and 4:1 
with results as shown in Figure 2¥*, 
By the combined use of the higher 
temperature and the higher acid to 
ore ratio, a more rapid extraction of 
manganese was obtained. After sev- 
eral hours, the percentage extraction 
of manganese was approximately the 
same in both experiments. The two 
experiments show that with an ade- 
quate amount of acid and with suit- 
able reaction conditions of time and 
temperature, open vessels can be used 
for leaching operations. 

Having produced manganous sul- 
phate liquors by single-stage leaching, 
it was necessary to determine if re- 
cycling would build up the concen- 
tration of manganese in the liquor. 
An aerated slurry of 200 gm 4.3% 
Mn ore and 10 gm pyrite in 800 ml 
2% v/v H,SO, (pH 0.4) was heated 
ina pachuca at 90°C. for four hours. 
After separating the leach liquor 
from the ore, some water was added 
to the liquor to compensate for that 
which adhered to the ore. The pH of 
the liquor was restored from approx- 
imately 1.0 to its original value of 0.4 
with about 25 gm conc. H,SO, and 
the liquor was recycled w ‘ith fresh 
ore and pyrite to the pachuca. The 
process was repeated for 11 cycles 
with the acidity of the leach liquor 
adjusted only for the first six cycles, 

Figure 3 shows that during the first 
Six cycles, when H,SO, additions 
were ‘made, there was a build up in 
the manganous sulphate concentra- 
tion of the leach liquor to about 100 
gm MnSO,/1. Although the average 
manganese extraction was 65%, in 
certain cycles the extraction was as 
high as 90%. Subsequent to the sixth 
cycle, the manganese concentration 
of the liquor increased slightly and 
then decreased. The decrease was 
probably due to a combination of the 
lower manganese extraction during 
the last few cycles and the dilution 
involved at each cycle in restoring 
the volume of the liquor to its initial 
value. During the cyclic leaching ex- 
periments, about two moles of sul- 
phuric acid were required for every 
mole of manganese that was leached. 

After six cycles, the leach liquor 
contained 3 gm Fe+++/1 and 7 gm 
Fe++/1l. An important fact noted 
during later cycles was that when 
the pH incre eased to above approxi- 
mately 1.5, the iron was Pye arene 
by hy droly sis; after cycle 8, no iron 
was detected in the leach liquor. 
Hence, the cyclic leaching of manga- 
nese from the low -gr rade pyrolusite 
ore produced a strong manganous 
sulphate solution that was free of iron. 
*In all graphs where pressures are not indicated, 


the experiments were done at atmospheric 
pressure in a pachuca. 
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Fyaporation of the recycled leach 
liquor gave a hydrated manganous 
a, The main impurities were 
Al, Ca, Mg, and Si, and these com- 
prise. d less than 4% by weight of the 
manganous sulphate. Manganous sul- 
phat monohydrate, the solubility of 
which decreases at high temperatures, 
was recovered by heating the re- 
cycled leach liquor to about 175°C. 
in a closed tube, using a thermal pre- 
cipitation method developed at the 
Mines Branch.* 


Study of the leaching reactions 


To elucidate the chemistry of the 
leaching process, a number of ex- 
periments were done with ore or pure 
MnO, using a glass pachuca and 
reaction temperature of 90°C. 


The effect of pyrite in the leaching 
process 


To demonstrate that the leaching 
of manganese involved more than a 
simple acid leach, the following ex- 
periment was done. An aerated slurry 
of ore in 2% v/v H,SO, was heated 
for two hours in the absence of pyrite. 
The extraction of manganese was 
about 7%. (Possibly this represents 
the percentage of the manganese that 
is present in the divalent state in the 
ore). As soon as pyrite was added, 
the percentage extraction of manga- 
nese from the pyrolusite ore was 
greatly increased as shown in Figure 
i 


Reactions of ferrous sulphate in 
the leaching process 

Since ferrous sulphate is probably 
the main leaching reagent in the pro- 
cess, a study was made of its relevant 
reactions and their comparative rates. 
The reaction by which manganese 
dioxide is dissolved by an_ acidic 
ferrous sulphate solution is usually 
written as (6) 


MnO, + 2FeSO, + 2H»SO, — MnSO, + 
Fes (SO,); + 2H:O. .(1) 


The rate of this reaction was very 
rapid at 90°C. as shown by Figure 5, 
with the dissolution of the manganese 
dioxide being complete in about one 
minute. In agreement with Equation 
1, the curves of Figure 5 show a 
change in pH and in ferrous ion con- 
centration during the dissolution of 
the manganese. 

Although the reaction of ferrous 
sulphate with manganese dioxide is 
fast, the formation of ferrous sulphate 


°In preparation for publication. 

°°When pyrrhotite, Fe,S., was substituted for 
pyrite, no significant change was observed in 
the leaching reaction rate. See also recent 
paper by Bjorling(11) on the autoclave 
leaching of manganese ore by pyrrhotite. 
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Figure 3—Cyclic extraction of manganese from ore in the presence of pyrite. 
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Figure 4—Effect of addition of pyrite on extraction of manganese from ore. 


The Canadian Journal of Chemical Engineering, February, 1958 39 





(GM/L) 


CONCENTRATION 


CONCENTRATION (GM/ L) 


Fe 


40 


5 gm MnOo 
45 gm Fe SQ, 7 H,0 
12 600 ml 2% H,SO, 
90°C 


wo 


Oo 2 4 ~ 6 8 10 
TIME (MINUTES) 


Figure 5—Rate of reduction of MnO: by ferrous sulphate. 
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Figure 6—Rate of oxidation of ferrous sulphate by air and by oxygen. 








from pyrite (5) is known to be slow. 


FeS, + 3140: + HO — FeSO, + H:S0, 

Ae aeons (2) 
For example, only 5% of the pyrite 
was oxidized in four hours when a 
slurry of minus 200 mesh pyrite in 
dilute H. SO, of pH 1.0 at 90°C. was 
aerated with oxygen. No significant 
change in the reaction rate was ob- 
served when an H,SO, solution of 
pH 0.4 was used. 

The ferrous sulphate formed by 
reaction 2 may be oxidized (+) by 
oxygen according to 
2FeSO, + 402 + HSO, — Fez (SO,); + 

tr 20. .(3) 
The removal of ferrous ion by this 
reaction might be expected to inter- 
fere with the leaching of pyrolusite 
by reaction 1. However, Figure 6 
shows that even when copper sulphate 
is present as a catalyst, the oxidation 
of ferrous sulphate by air or Oxy gen 
is relatively slow. Hence the air oxi- 
dation of ferrous ion in acidic solu- 
tion is not likely to impede the leach- 
ing operations. 

Although equations 2 and 3 indi- 
cate that ferric sulphate could be 
formed when pyrite is oxidized, no 
ferric ion was detected when acidic 
slurries of pyrite at 90°C. were 
aerated at atmospheric pressure. Any 
ferric ion that was formed must have 
been reduced again by the residual 
pyrite, possibly by a reaction such 
as (2) 

FeS.+ 7Fe2 (SO,4)3 +'8H20 — 

15FeSO, + 8H2SO,. .(4) 
The rate of reduction of a ferric 
sulphate solution by minus 200 mesh 
by rite is shown g! raphically by Figure 

. Although this is a slow reaction, it 
sight be, for liquors containing ferric 
sulphate, the main source of ferrous 
ion for leaching pyrolusite ore. 


The effect of acid concentration 
and of ferric sulphate addition on 
the leaching process 


In the first experiment, a slurry of 
MnO, and pyrite in an H,SO, solution 
of pH 0.4 at 90°C. was aer vined with 
oxygen. Figure 8 shows that the man- 
ganese was dissolved at a uniform 
rate; after four hours approximately 
80% of the manganese had reacted 
with the py rite to form manganous 
and ferric ions. The initial rise in the 
pH of the solution was probably due 
to the impurity siderite, FeCO,, that 
was present in the pyrite and the 
subsequent progressive rise in the pH 
was likely due to the dissolution of 
the manganese dioxide. 

When dilute H,SO, of pH 1.0 was 
substituted for acid of pH 0.4 in the 
above experiment, the results shown 


The Canadian Journal of Chemical Engineering, February, 1958 


——— OO OOOO ee eee 
A A 





in | 
dili 
of | 
ren 
fer! 
sho 
pro 
tio! 


to 

anc 
of 

inc 
Mr 
wa 
fer! 
wit 
tin! 
aft 


Au 


sul 


by 
tiv! 
ob 
of 

aut 
an 

H« 
vic 
lea 
atu 
sul 


py 
in 

11 
of 
de 
i01 
of 
py 
ne 
be 
30 
sul 
pr 
pr 
ob 
sal 
OX 
an 


of 
qu 
ro 


iv 





Ne 


hd 


————— a NEE EEE rR ea SSS SSS 


in Figure 9 were obtained. With the 
dilute acid of this experiment, most 
of the acid was soon used up and the 
remaining was insufficient to keep 
ferric ion in solution. The curves 
show that for pH values above ap- 
proximately 1.7, the rate of dissolu- 
tion of manganese was very slow. 
The effect of adding ferric sulphate 
to an acidic aerated slurry of MnO, 
and pyrite, as shown by a comparison 
of Figure 8 with Figure 10, was to 
increase the rate of dissolution of the 
MnO,. Evidently the ferric sulphate 
was reduced by the pyrite to form 
ferrous sulphate which reacted rapidly 
with the MnO,. The ferric ion con- 
tinued to be reduced by the pyrite 
after all the MnO, had dissolved. 


Autoclave-produced acidic ferrous 
sulphate solutions 


Since the leaching of pyrolusite ore 
by pyrite consumes acid, an alterna- 
tive procedure was investigated. The 
objective was to adapt the oxidation 
of an aqueous slurry of pyrite in an 
autoclave (7-19) for the production of 
an acidic ferrous sulphate solution. 
Hoak and Coull “') have shown pre- 
viously that pyrolusite ore can be 
leached rapidly at ambient temper- 
atures with acidic solutions of ferrous 
sulphate. 

In one experiment, a slurry of 
pyrite in 2% v/v H,SO, was heated 
in an autoclave at a temperature of 
110°C. and under an oxygen pressure 
of 300 psig. The acid was added to 
decrease the precipitation of ferric 
ion by hydrolysis in the initial stages 
of the reaction. About 60% of the 
pyrite was oxidized in five hours with 
nearly 70% of the pyritic sulphur 
being changed to sulphate sulphur and 
30% being recovered as elemental 
sulphur. This autoclave treatment 
produced an acidic liquor containing 
predominantly ferric sulphate. To 
obtain a ferrous sulphate liquor, the 
same experiment was repeated but the 
oxygen was released after three hours 
and replaced by nitrogen. 

Figure 11 shows that the treatment 
of pyrite in an autoclave produced 
quickly a liquor containing both fer- 
rous and ferric sulphate. When a high 
oxygen pressure was present, the iron 
was changed mainly to the ferric state. 
When the oxygen was replaced by 
nitrogen, the reauction of the ferric 
ion by pyrite produced an_ acidic 
liquor containing approximately 200 
gm FeSO,/, 60 gm Fe,(SO,),/1 and 

70 gm H,SO, /\. This oe 
liquor is similar in composition to 
pickle liquor which has been reported 
‘) to be suitable for leaching pyro- 
lusite ore at room temperature. The 
autoclave process appears to be a con- 
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Figure 8—Rate of dissolution of MnO, in 2% v/v H.SO, in the presence of pyrite. 
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pyrite and ferric sulphate. 


venient method of producing an 
acidic ferrous sulphate liquor for 
leaching purposes. 


Summary and Conclusions 


Manganese was extracted from a 
pyrolusite ore (4% Mn) by heating 
an aerated aqueous slurry of ore and 
pyrite. A rapid extraction of manga- 
nese was obtained by leaching the ore 
in an autoclave at high oxygen pres- 
sures and high temperatures. It was 
found that satisfactory extraction 
could be obtained, at a slower rate, 
by using open vessels and relatively 
low temperatures. A temperature of 
90°C. was used for most of the 
experimental work since adequate 
extractions could be obtained in a 
relatively short time, e.g. manganese 
extractions of up to 90% were ob- 
tained within four hours. For com- 
mercial operations, prolonged leach- 
ing periods at lower temperatures 
might be more economical. 

Recycling the leach liquor at con- 
trolled acidity to leach additional ore 
gave a concentrated solution contain- 
ing about 100 gm MnSQ,/1. Approxi- 
mately two moles of sulphuric acid 
were required to dissolve one mole of 
pyrolusite from the ore. Evaporation 
of the recycled solution gave an iron- 
free hy drated manganous sulphate of 
which the main impurities, Al, Ca, 
Mg and Si, comprised less than 4% 
by weight of the product. Manganous 
sulphate was recovered using substan- 
tially less heat when evaporation was 
replaced by heating the recycled 
solution to about 175°C. to precipi- 
tate manganous sulphate monohy- 
drate. 


Although many chemical reactions 
are undoubtedly involved, the leach- 
ing of pyrolusite ore by this process 
is essentially the reduction of MnO, 
by ferrous ion to acid- solution diva- 
lent manganese. The ferrous ion is 
obtained from the reaction of pyrite 
with air or with ferric ion. During 
the dissolution of manganese, the iron 
is alternately oxidized by the pyro- 
lusite and reduced by the pyrite, as 
shown by 


MnO, 
Fett === Fet* 
FeS, 
i.e. the ferrous ion is continuously 
regenerated during the leaching pro- 
cess. 

By heating an aqueous slurry of 
pyrite at 110°C. in an autoclave, first 
under oxygen pressure and then in 
the absence of oxygen, a concentrated 
acidic solution of ferrous sulphate 
(200 gm FeSO,/1) was obtained. 
Hoak and Coull “) found that solu- 


tions of similar composition were 
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suitable for leaching pyrolusite ore at 
room temperature. Hence, for large- 
scale leaching operations with the 
low-grade ore used in this investiga- 
tion, two leaching procedures might 
be feasible. In the one process, the 
leaching would be done by heating 
aerated slurries of ore, pyrite, acid 
and water in pachucas. In the altern- 
ative process, pyrite and water, in 
amounts small in comparison with 
that of the ore, would be treated in 
autoclaves to produce acidic solutions 
of ferrous sulphate for the subsequent 
leaching of ore at ambient temper- 
atures. 
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